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ABSTRACT 

This  paper  presents  the  development  of  a  graphic  computer  simulation  of  a 
Jackson  network.  In  particular,  the  simulated  network  is  a  queueing  system  consisting 
of  two  servers  in  tandem,  one  customer  entrance  point,  two  potential  customer  exit 
points,  and  one  potential  customer  feedback  path. 

This  system  is  modeled  as  a  discrete-event  simulation  implemented  on  an  IBM 
Personal  Computer  (IBM  PC).  The  IBM  PC  requires  a  color  monitor,  a  color/graphics 
adapter  card,  at  least  one  disk  drive,  and  a  minimum  of  128  kilo-bytes  of  random 
access  memory  (RAM). 
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I.  INTRODUCTION 

This  paper  presents  the  development  of  a  computer  simulation  of  a  specific 
example  of  a  Jackson  network. 

In  a  typical  computer  simulation,  a  series  of  mathematical  and  logical 
relationships  are  used  to  create  a  model  of  some  real  world  system.  The  model  is  then 
exercised  over  some  period  of  time,  and  data  are  collected  to  estimate  the  desired 
characteristics  of  the  real  word  system.  When  the  simulation  run  is  completed,  the 
desired  statistics  are  printed  out. 

What  makes  this  simulation  different,  is  that  it  also  provides  the  user  with  a 
visual  picture  of  what  is  occuring  in  the  queue  network  as  the  simulation  progresses.  It 
accomplishes  this  by  actually  displaying  on  the  computer  screen  the  arriving  and 
departing  customers  and  showing  how  they  move  through  the  network.  This  portrayal 
is  possible  because  the  simulation  gives  the  observer  the  impression  of  real-time 
passage. 

Because  Jackson  networks,  and  their  associated  statistical  quantities  of  interest 
can  be  solved  in  closed  mathematical  form,  there  is  little  gained  in  this  respect  from 
simulation.  However,  by  allowing  a  student  to  actually  observe  the  inner  workings  of 
the  queue  network  as  it  progresses  in  time,  insights  into  the  process  can  be  obtained 
that  are  not  possible  from  conventional  instructional  tools  alone. 

This  paper  is  the  third  known  effort  at  the  Naval  Postgraduate  School  to  develop 
visual  simulations  of  various  stochastic  processes  utilizing  the  personal  computer. 
Prior  undertakings  are  a  simulation  of  various  models  that  are  based  on  an  underlying 
Poisson  process  [Ref  1],  and  a  simulation  of  a  Machine- Repairman  Model  [Ref  2]. 
These  two  papers  provided  many  ideas  for  the  writing  of  the  Jackson  network 
simulation  and  served  as  an  outline  for  this  paper. 

The  remainder  of  this  paper  describes  the  Jackson  network  model  and  its 
implementation  on  the  IBM  Personal  Computer  (IBM  PC).  Chapter  IV  is  a  user's 
guide  which  contains  instructions  for  actually  running  the  simulation.  A  program 
listing   is  contained  in  Appendix  A. 


II.  THE  JACKSON  NETWORK  MODEL 

A.       GENERAL  DISCUSSION 

In  a  network  of  queues,  a  customer  may  receive  service  from  one  or  more 
service  centers,  possibly  revisiting  some  or  all  of  the  service  centers  many  times  before 
finally  exiting  the  network.  With  no  additional  assumptions  made,  this  type  of  queue 
network  can  be  extremely  difficult,  if  not  impossible  to  analyze  in  closed  mathematical 
form.  In  these  cases,  computer  simulation  may  be  the  only  alternative  [Ref  3:  p.  225]. 

One  type  of  network  that  does  permit  analysis  without  the  need  for  simulation  is 
the  Jackson  network.  A  Jackson  network  is  a  collection  of  queues  that  possess  the 
following  characteristics  [Ref  4:  p.  456]: 

•  The  network  consists  of  N  service  centers  numbered  1,2,...,N.  Service  center  i 
contains  c-  identical  servers. 

•  Customers  from  outside  the  network  arrive  at  service  center  i  according  to  a 
Poisson  process  with  rate  A-,  i  =   1,2,. ..,N. 

•  After  receiving  service  at  center  i,  a  customer  leaves  the  network  with  probability 
PjQ  >  0  or  goes  instantaneously  to  service  center  j  with  probability  ?■:,  ^P-  =  1. 
These  probabilities  are  independent  of  the  history  of  the  network. 

•  Customers  arriving  at  service  center  i  (from  inside  or  outside  the  network)  are 
served  according  to  a  First  Come  First  Serve  (FIFO)  queue  discipline.  The 
service  times  are  independent,  identically  distributed  exponential  random 
variables  with  mean  i/ji-,  i  =   1,2,.. .,N. 

•  There  is  unlimited  waiting  space  at  every  service  center. 

•  The  effective  arrival  rate  at  each  service  center  taking  into  account  arrivals  from 
both  outside  and  inside  the  network  must  be  less  than  the  corresponding  service 
rates. 

Two  further  characteristics  of  a  Jackson  network  result  from  the  requirement  that 

customers  from  outside  the  network  arrive  at  the  service  centers  according  to  a  Poisson 

process.    This  implies  that  the  arrival  process  must  have  independent  and  stationary 

increments.     An  arrival   process   is   said   to   possess   independent   increments  if,   for 

example,  the  number  of  arrivals  that  occur  between  times  5  and  10  are  independent  of 

the  number  of  arrivals  between  times  35  and  40.    An  arrival  process  has  stationary 

increments  if  the  distribution  of  the  number  of  events  which  occur  in  any  interval  of 

time  depends  only  on  the  length  of  the  time  interval.    This  is  the  same  as  saying  that 


^Although  this  paper  will  refer  to  customers,  implying  people,  the  flow  through  a 
queue  network  could  just  as  easily  be  an  object  or  a  job. 


there  is  no  lime  of  day  cflect  on  the  arrival  rate  of  customers.  One  further  assumption 
of  a  Jackson  network  is  that  the  servers  at  the  difTerent  service  centers  work 
independently  and  do  not  assist  one  another,  even  if  one  is  idle. 

B.       THE  MODELED  JACKSON  NETWORK 

The  simulation  described  in  this  paper  is  a  specific  example  of  a  Jackson  network 
where  N  =  2,  c^  =  1.  c^  =  1,  l^  >  0,  X^  =  0,  \i^  >  0.  ^2  >  0,  P,,  =  0,  \\^  =  0, 
P  >  0,  P  >  0,  Pj2  =  1  -  ^10'  ^^^  ^21^  ^  '  ^20-  ^'^  pictorial  representation  of 
this  network  is  provided  in  Figure  2.1.  In  words,  a  customer  arriving  from  outside  the 
network  must  first  go  to  service  center  1  where  he  enters  a  FIFO  queue.  Upon 
completion  of  service  at  service  center  1,  the  customer  can  either  exit  the  system  with 
probability  P^^  or  proceed  to  ser\'ice  center  2  with  probability  1  -  P^^,  where  he  again 
enters  a  FIFO  queue.  Upon  completion  of  service  at  service  center  2,  the  customer 
can  either  exit  the  network  with  probability  P^^  or  loop  back  to  service  center  1  with 
probability  1  -  P^q  again  entering  the  FIFO  queue.  In  this  configuration  a  customer 
can  loop  from  service  center  1  to  service  center  2  and  back  to  service  center  1 
numerous  times  before  exiting  the  network.    Each  service  center  has  only  one  server. 
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Figure  2.1     Jackson  Network  With  Two  Service  Centers. 
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It  is  arguable  whether  any  real  world  system  would  ever  meet  all  of  the 
requirements  of  a  Jackson  network.  But  then  that  is  always  the  predicament  faced  by 
the  modeler;  making  the  model  represent  the  real  world  without  making  it 
unreasonably  complex.  There  are  numerous  real  world  systems  that  might  be 
reasonably  modeled  by  this  simple  Jackson  network.  For  example,  the  network  might 
describe  the  flow  of  a  piece  of  equipment  through  a  repair  facility.  Upon  entering 
service  center  1,  the  repairman  evaluates  the  piece  of  equipment  and  attempts  to  repair 
it.  If  it  can  not  be  repaired,  it  leaves  the  network.  If  it  is  repairable,  it  is  repaired  and 
proceeds  to  the  Quality  Assurance  Shop  (service  center  2).  At  Quality  Assurance,  the 
repaired  item  undergoes  testing  and  is  either  classified  as  fit  or  unfit.  If  it  is  determined 
fit,  it  departs  the  network  and  returns  to  the  field.  If  it  is  determined  unfit,  it  is 
returned  to  service  center  I  for  additional  repair. 

Another  possible  example  would  be  a  computer  system.  Jobs  arrive  at  the 
central  processing  unit  (service  center  I)  of  a  computer.  When  a  job  finishes  CPU 
service  it  either  departs  the  computer  system  or  requires  Input/Output  service  (service 
center  2).  When  a  job  completes  lO  service  it  either  leaves  the  computer  system  or 
returns  to  the  CPU  for  additional  service. 

For  this  Jackson  network  model,  the  state  of  the  system,  N(t)  =  n,  is  defined  as 
the  total  number  of  customers  in  the  system  both  waiting  and  being  served  at  time  t. 
This  state  variable  is  further  broken  down  into  the  total  number  of  customers  waiting 
and  being  served  by  service  center  I,  N.(t)  =  n  ,  and  the  total  number  of  customers 
waiting  and  being  served  by  service  center  2,  'N^{i)  =  n^.  It  follows  that 
N(t)  =  N-(t)  +  N2(t).  Because  arrivals  and  departures  are  discrete  events,  the  value 
of  N(t)  will  always  be  zero  or  a  positive  integer.  However,  because  there  is  unlimited 
waiting  space  in  front  of  the  service  centers,  there  is  no  upper  limit  on  the  value  that 
N(t)  can  assume.   That  is,  n  can  be  any  integer  value  between  0  and  positive  infinity. 

C.       ANALYSIS  OF  JACKSON  NETWORKS 

It  is  not  the  intent  of  this  section  to  develop  in  detail  the  theoretical  analysis  of 
Jackson  networks.  However,  a  general  discussion  will  be  included  for  completeness. 
The  interested  reader  should  refer  to  the  reference  list  and  the  bibliography  for  books 
on  queueing  theory,  paying  particular  attention  to  the  chapters  on  networks  of  queues. 
Of  particular  interest  might  be  the  paper  written  by  J.R.  Jackson  for  whom  this  type  o[ 
network  is  named. 
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By  imposing  the  restrictions  described  in  Section  I-A  on  a  network,  Jackson 
showed  that  in  the  steady  state  each  service  center  (say  center  i)  in  the  network  behaved 
as  if  it  were  an  independent  M/M/c  queue  with  service  rate  u-  and  Poisson  input  rate 
Rj.  The  input  rate,  R-,  can  be  considered  the  eflective  arrival  rate  at  service  center  i. 
The  effective  arrival  rate  to  a  service  center  is  the  sum  of  the  Poisson  arrivals  from 
outside  the  system  (exogenous  customers),  plus  the  arrivals  from  any  other  service 
center  in  the  network.  The  effective  arrival  rates  can  be  found  by  solving  the  following 
set  of  simultaneous  equations: 

^[=\  +  ^RjPj-     i  =  I,2,...,N.  (eqn  2.1) 

In  order  to  prevent  an  unstable  or  explosive  queue  from  occurring  at  any  of  the 
service  centers,  it  is  required  that  R-  <  c-u-.  In  words,  the  effective  arrival  rate  of 
customers  at  service  center  i  must  be  less  than  the  cumulative  service  rate  of  the  servers 
at  service  center  i. 

The  implications  of  Jackson's  findings  are  that  a  complex  network  can  be 
decomposed  into  a  number  of  simpler  M/M/c  sub-systems.  These  sub-systems  can 
then  be  analyzed  using  Little's  formula  and  the  standard  results  for  the  M/M/c  queue. 

An  interesting  aspect  of  this  result  is  that  it  can  be  proved  that  allowing  feedback 
in  the  network  destroys  the  Poisson  flow  of  arrivals  [Ref  4:  p.  460].  However 
Jackson's  results  show  that  even  for  networks  that  allow  feedback,  the  service  centers 
behave  as  if  they  were  fed  totally  by  Poisson  arrivals,  w^hen  in  fact  they  are  not. 
Statisticians  have  been  unable  to  provide  an  explanation  as  to  why  this  is  true  [Ref  5: 
p.  330]. 

D.       MODEL  ESTIMATIONS 

The  purpose  of  any  simulation  is  to  gain  some  insight  into  the  real  world  system 
being  modeled.  This  insight  is  usually  gained  by  estimating  certain  quantities 
pertaining  to  the  systems  performance.  As  a  part  of  the  graphical  output,  the 
simulation  provides  a  table  of  estimates  for  a  number  of  these  quantities.  Because 
Jackson  networks  can  be  analyzed  in  closed  form,  the  graphical  output  will  also  include 
the  theoretical  values  for  these  same  quantities.  This  will  provide  the  user  with  a 
measure  of  how  well  the  computer  simulation  is  modeling  the  real  world  system. 
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The  following  equilibrium  performance  quantities  will  be  both  estimated  and 
determined  theoretically  by  the  simulation: 

R-  =  the  effective  arrival  rate  of  customers  at  service  center  i  (i  =   1,2). 

S;  =  the  expected  number  of  times  a  customer  will  pass  through  service 
center  i  before  exiting  the  network  (i  =   1,2). 

L;  =  the  expected  number  of  customers  both  waiting  and  being 
served  at  service  center  i  (i=   1,2). 

L     =  the  expected  number  of  customers  in  the  network,  both  waiting 
and  being  served. 

W|  =  the  expected  amount  of  time  a  customer  spends  waiting  and 
being  served  at  service  center  i  (i  =   1,2). 

W     =  the  total  expected  amount  of  time  a  customer  spends  in  the  network,  both 
waiting  and  being  served. 

Z     =  the  proportion  of  time  service  center  one  is  busy  and  service 
center  two  is  idle. 

The  theoretical   values  will  be   represented  by  the  appropriate   capital  letter.     The 

estimated  values  will  be  represented  by  the  appropriate  capital  letter  with  a  hat  over  it. 

The  following  quantities  are  parameters  of  the  model  which  can  be  selected  by 

the  user  prior  to  running  the  simulation  and  are  used  in  determining  the  theoretical 

quantities  just  listed.    The  parameter  names  are  consistent  with  those  used  in  the 

simulation  model. 

X  =  the  Poisson  arrival  rate  of  exogenous  customers  arriving  at  service  center  1 
(previously  referred  to  as  X^). 

]i^  =  the  exponential  service  rate  for  service  center  (i=  1,2). 

P;  =  the  probability  of  exiting  the  network  after  completing  service  at  center  i, 
i=  1,2  (previously  referred  to  as  P-  ). 

1  -  P.  =  the  probability  of  proceeding  to  service  center  2  after  completing 
service  at  service  center  1  (previously  referred  to  as  P^^^)- 

1  -  P„  =  the  probability  of  proceeding  to  service  center  1  after  completing 
service  at  service  center  2  (previously  referred  to  as  P^^). 

For  use  in  estimating  some  of  the  performance  quantities,  X  and  ^-  will  also  be 
estimated  by  the  simulation  model. 
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The  values  of  R.  and  R^,  the  theoretical  efTective  arrival  rates  at  service  center  1 
and  service  center  2,  are  found  by  solving  equation  2.1,  namely: 

R^  =  X  +  RgCl-Pg)  (eqn2.2) 

R^  =  0  +  Ri(l-P,),  (eqn2.3) 

which  results  in: 

R^  =  Vil-d-PiXl-P^)}  (eqn2.4) 

R^  =  (l-Pi)Ri.  (eqn2.5) 

The  estimated  value  of  R^  is  computed  by  recording  the  total  number  of 
customer  arrivals  at  service  center  1,  both  exogenous  customers  and  feedback 
customers  from  service  center  2,  and  dividing  by  the  total  elapsed  time  that  the 
simulation  has  run.  The  estimated  value  of  R^  is  calculated  in  a  similar  manner  except 
that  there  is  only  one  source  of  arrivals  to  service  center  2,  namely  those  customers 
departing  service  center  1  who  do  not  exit  the  network. 

The  estimated  value  of  X,  the  Poisson  arrival  rate  of  exogenous  customers 
arriving  at  service  center  I  is  computed  by  recording  the  number  of  exogenous 
customers  arriving  at  service  center  1  and  dividing  by  the  total  elapsed  time  that  the 
simulation  has  run. 

The  estimated  value  of  ]i  ,  the  exponential  service  rate  for  service  center  1  is 
computed  by  recording  the  total  number  of  customers  completing  service  at  service 
center  1,  and  dividing  by  the  total  amount  of  time  service  center  1  was  busy  (  NOT  the 
total  elapsed  time  of  the  simulation).  The  estimated  value  of  H-  is  computed  in  a 
similar  fashion. 

The  value  of  L^^  and  L^,  the  theoretical  expected  number  of  customers  either 
waiting  or  being  served  at  service  center  1  and  2,  are  found  by  utilizing  Little's  law  and 
the  results  for  a  M/M/1  queue: 

Lj  =  R,/{M,-R,)  (eqn  2.6) 
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L^  =  R^/CM^-R^).  (eqn  2.7) 

It  follows  that  L,  the  expected  number  of  total  customers  in  the  network,  both 
waiting  and  being  served  is  simply  the  sum  of  L^^  and  L  : 

L  =  L^  +  L^.  (eqn  2.8) 

/X       A  A  A         A        A  A 

The  estimates  L^  .L^  and  L  are  found  by  substituting  R  ,  R  ,  ^  ,  and  ^  into  the 
above  equations  to  obtain: 

L^  =  R,/(?,-R,)  (eqn  2.9) 

L^  =  Mlle-Rg)  (eqn  2.10) 

A  A  A 

L  =  L^  +  Lg.  (eqn  2.11) 

The  values  of  W^^  and  W  ,  the  theoretical  expected  time  a  customer  spends 
waiting  and  being  served  at  service  center  1  and  2,  are  also  found  by  utilizing  Little's 
law  and  the  results  for  a  M/M/1  queue: 

W^  =  L^/  R^  (eqn  2.12) 

W^  =  Lg/  Rg.  (eqn  2.13) 

The  value  of  W,  the  total  expected  time  a  customer  spends  in  the  network,  both 
waiting  and  being  served  is  also  computed  using  Little's  formula: 

W  =  L/X.  (eqn  2.14) 

The  estimates  W^^,  W^,  and  W  are  found  by  substituting  L,  L^,  L^,  R^,  R^,  and  X 
into  the  above  equations  to  obtain: 

/N  A        A 

W^  =  L^/Ri  (eqn  2.15) 
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/\  A        A 

W,  =  LJK  (eqn2.16) 


A  A   A 

W  =  L;X.  (eqn2.17) 


The  value  of  Z,  the  theoretical  long  run  proportion  of  time  that  service  center  1 
is  busy  and  service  center  2  is  idle,  is  found  by  using  a  result  commonly  referred  to  as 
Jackson's  Theorem.  This  theorem,  when  simplified  for  the  situation  where  each  service 
center  has  exactly  one  server,  states  that  for  a  Jackson  network  with  Rj  <  u-  and  with  k 
service  centers  the  limiting  probability  is  given  by: 


Prob(n-  at  service  center  i,  i=  1,2,. ..k)  =    7t(nj^,n2,...nj^)  =  (eqn  2.18) 

k 

n 

i  =  l 


n{n^}n{n^)...n{n^)  =  n{R-J}i-f-^{l-R-l]i-). 


Where  n-  is  the  number  of  customers  waiting  and  being  served  at  service  center  i,  and 
7i(n-)  is  the  long  run  probability  of  finding  n-  customers  at  service  center  i.  The 
theorem  says  that  the  limiting  probability  is  equal  to  the  product  of  the  marginal 
probabilities  for  the  individual  M/M/1  queues  comprising  the  network. 

To  determine  the  value  of  Z,  it  is  necessary  to  determine  the  probability  that 
there  are  no  customers  at  service  center  2  and  one  or  more  at  service  center  1.  This 
probability  is  determined  as  follows: 

Z  =  (1 -7r(n^  =  0)}7r(n2  =  0)  (eqn  2.19) 

=  {l-(Ri/Hi)0(l-Ri/^i)}{(R2/H2)^(l-R2/H2)}  (eqn  2.20) 

=  (R^/^^)(l  -  R^/Mp.  (eqn  2.21) 

A  A       A      A  A 

The  estimate  Z  is  found  by  substituting  K^,  R^,  H^,  and  ^^  into  equation  2.21  to 

obtain: 

A  A     /\  A      A 

Z  =(R^/M^)(l-R2/iip.  (eqn  2.22) 
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The  value  of  S^,  the  theoretical  expected  number  of  times  a  customer  will  pass 
through  service  center  1  before  exiting  the  network  is  found  by  recognizing  that  the 
underlying  distribution  is  geometric.  For  service  center  1,  it  can  be  seen  that  the 
following  relationship  holds,  where  N^^  equals  the  number  of  times  a  customer  visits 
service  center  1  before  exiting: 

P(N^=1)=  {P^  +  (l-Pi)?^}  (eqn2.23) 

P(N^=2)  =  {P^  +  (1-P^)P2}((1-P^){I-PP}1  (eqn  2.24) 

P(N^=3)  =  {P^  +  (1-P^)P2}{(I-P^)(1-P2)}2  (eqn  2.25) 


P(N^=n)  =  {P^+  (1-P^)P2}{(1-P,)(1-P2)}"-'   ,  n=  1,2,3....  (eqn  2.26) 

Expanding  the  right  hand  side  of  equation  2.26, 

P(N  =  n^)  =  (P^+P2-P,P/(l-P,-P2  +  PiP2)^-'  (eqn  2.27) 

Letting  9  =   1  -  Pj^  -Pg  +  PiP?'  ^^  follows  that  the  quantity  raised  to  the  first  power  in 
equation  2.27  is  equal  to  1  -  0.    Equation  2.26  can  be  rewritten  in  the  form: 

P(N  =  n)  =  (i-0)(e)^-\   n  =  1,2,3....  (eqn  2.28) 

Which  is  recognized  as  a  geometric  distribution  with  an  expected  value  equal  to: 

Si  =  E(N^)  =   1/(P^  +  P^  -  P.Pp.  (eqn  2.29) 
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A  similar  type  of  analysis  can  be  done  to  find  the  theoretical  value  of  S^.  Letting 
N  equal  the  number  of  times  a  customer  visits  service  center  2  before  exiting  and 
0  =    (1-P  Xl-P^)  as  before,  it  can  be  seen  that  the  following  relationships  hold: 

P(N2  =  0)  =  P^  (eqn2.30) 

P(N2=  1)  =  ePj  +  (l-Pi)P2  (e^n  2.31) 

P(N2=2)  =  02p^  +  Q\l-?^)?^  (eqn  2.32) 


P(N2=n)  =  e^P^.  +  e^-Hl-P^)?^.  (eqn  2.33) 

Taking  the  expected  value: 

E(N2)  =  OPj  +  le^P^  +2e2p^  +  ...+  ne^P^  (eqn  2.34) 

+  l(l-Pi)P2  +  2eHl-Pj)P2  +  3G2(1.Pj)P2  +  ...+  n0^-^(l-P|)P2 

=  P^0{1  +  20^+  302  +....  +  n0"-i}  (eqn  2.35) 

+  (1-P^)(P2){1  +  20^  +  302  +....  + n0"-^}. 

Noting  that  the  geometric  sequences  converge  in  the  limit  to  (1/1-0)2,  equation  2.35 
can  be  rewritten  as: 

E{N^)  =  P,0(l/l-0)2  +(l-P^)(P2)(l/l-0)2.  (eqn  2.36) 

The  final  result  after  manipulation  of  equation  2.36  is: 

S^  =  E{y^)  =  il-?^)/{?^+P^.?^P^)  =  (1-P,)S^.  (eqn  2.37) 
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The  estimates,  S^^  and  S^,  are  found  by  using  the  definition  of  expected  value  for 
a  discrete  random  variable.  Let  Pj(n)  be  a  shorter  notation  for  the  probability  that  a 
customer  is  served  exactly  n  times  at  service  center  i  before  he  departs  the  network. 
The  estimate  P;(n),  is  determined  by  recording  how  many  customers  were  served 
exactly  n  times  by  service  center  i  before  departing  the  network  and  then  dividing  this 
number  by  the  total  number  of  customers  who  have  departed  the  network.  The 
estimates  of  S    and  S^  are  then  determined  by: 

A  00    A 

Sj  =  SnP^(n)  (eqn  2.38) 

n=0 

A         29  A 

S^  =  InP^Cn).  (eqn  2.39) 

n=0 


It  should  be  noted  that  the  simulation  program,  JACKQUE,  is  incapable  of 
explicitly  tracking  an  individual  customer  for  more  than  nine  visits  to  either  service 
center.  All  customers  who  exit  the  network  having  visited  either  service  center,  ten  or 
more  times,  are  counted  as  if  they  had  visited  that  center  exactly  ten  times.  This 
results  in  S^^  and  S^  being  under-estimates  of  the  theoretical  values.  Section  III  will 
discuss  the  implications  of  this  model  limitation. 

Whenever  a  quantity  is  estimated  through  simulation,  the  estimate  will  be  a  point 
estimate  of  the  true  value.  This  point  estimate,  being  a  random  variable,  will  have 
some  variation  associated  with  it.  As  the  simulation  proceeds  in  time,  the  point 
estimate  will  fluctuate  around  the  true  value.  At  any  particular  point  in  time  the  point 
estimate  may  or  may  not  be  close  to  the  true  value.  To  overcome  this,  a  simulation 
model  is  typically  run  a  large  number  of  times  using  some  stopping  rule.  The  ending 
estimates  for  each  run  are  then  averaged  to  get  more  accurate  estimates.  Some  sort  of 
variance  reduction  technique  like  Jack-Knifing  or  Boot-Strapping  is  also  typically 
employed. 

In  JACKQUE  however,  the  simulation  model  is  exercised  only  once,  then 
terminated.  For  this  reason,  at  any  particular  point  in  time,  the  estimates  might  not 
coincide  closely  with  the  theoretical  values.  This  is  particularly  true  of  the  estimates 
for  L,  W,  and  Z  which  are  algebraic  combinations  of  other  point  estimates.    For 

A     >S 

example,  the  quantity  W^^  is  estimated  using  the  quotient  L^/R^.  If  at  the  point  m  time 
the  estimate  is  made,  L^  is  larger  and  R^  is  smaller  than  their  true  values,  the  resulting 
estimate  of  Z  will  be  very  much  larger  than  it's  true  value. 
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Since  the  primary  purpose  of  JACKQUE  is  not  estimation,  this  is  not  a  serious 
shortcoming  of  the  model. 
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III.  THE  SIMULATION 

A.       GENERAL 

The  Jackson  network  is  modeled  as  a  discrete-event  simulation  model.  The 
simulation  is  both  dynamic  and  stochastic  in  nature.  A  discrete  system  (vice 
continuous)  is  one  in  which  the  state  variables  can  only  change  at  a  countable  or  finite 
number  of  points  in  time.  These  points  in  time  are  the  ones  at  which  an  event  occurs. 
For  the  Jackson  network,  the  state  of  the  system  N(t),  can  only  change  as  a  result  of 
an  exogenous  arrival,  or  a  service  completion  by  either  service  center. 

A  dynamic  simulation  (vice  static)  is  one  that  changes  or  evolves  in  time.  The 
Jackson  network,  which  is  driven  by  inter-arrival  times  and  service  times,  is  obviously 
of  this  nature. 

A  stochastic  simulation  model  (vice  deterministic)  is  one  that  contains  one  or 
more  random  variables.  In  other  words,  some  characteristic  of  the  system  is  random  in 
nature.  For  the  Jackson  network,  both  the  exogenous  customer  arrival  rate,  and  the 
service  rates  are  random  variables  generated  from  an  underlying  exponential 
probability  distribution. 

A  dynamic  simulation  requires  some  mechanism  for  keeping  track  of  the  passage 
of  elapsed  time,  allowing  the  program  to  know  when,  and  in  what  order  events  will 
occur.  In  the  case  of  the  Jackson  network,  we  need  to  keep  track  of  when  arrivals  and 
service  completions  occur.  The  variable  in  the  simulation  model  that  performs  this 
task  is  called  the  simulation  clock.  The  unit  of  time  for  the  clock  is  usually  not 
explicitly  stated,  but  needs  to  be  consistent  with  the  units  of  the  input  parameters. 

There  are  two  basic  approaches  for  advancing  the  simulation  clock.  It  can  be 
accomplished  using  a  next-event  time  advance,  or  a  fixed  increment  time  advance.  In 
both  cases  the  clock  is  initialized  to  zero  and  the  times  of  all  future  events  are 
determined. 

With  a  fixed-increment  clock,  the  clock  is  advanced  in  increments  of  At.  After 
each  At  increment,  a  check  is  made  to  determine  if  any  events  should  have  occured.  If 
one  or  more  events  should  have  occured,  these  events  are  considered  to  have  occurred 
at  the  end  of  the  interval,  and  the  system  state,  N(t),  is  updated  accordingly.  Two 
obvious  problems  with  a  fixed-increment  clock  are  that  event  times  appear  to  be  larger 
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than  they  should  be  due  to  the  rounding  up  eflect,  and  if  more  than  one  event  occurs 
in  At  it  might  not  be  possible  to  decide  in  which  order  the  events  occurred. 

With  a  next-event  clock,  the  successive  jumps  of  the  clock  are  unequal  in  size, 
and  correspond  to  the  amount  of  time  till  the  next  event  occurs.  After  the  simulation 
clock  is  advanced  to  the  first  of  these  future  events,  the  state  of  the  system  N(t)  is 
updated  to  account  for  the  occurrence  of  the  event,  and  times  of  future  events  are 
generated  as  necessary.  The  next-event  clock  is  used  by  most  simulation  languages 
because  it  does  not  introduce  bias  into  the  estimates,  and  by  skipping  over  periods  of 
inactivity  the  program  runs  faster. 

For  purposes  of  the  graphical  display,  JACKQUE  uses  di  fixed-increment  clock. 
The  clock  is  increased  in  increments  of  At  =  I.  This  gives  the  viewer  of  the  graphics 
the  illusion  of  real  time  passage.  This  type  of  clock  is  also  necessary  so  that  the  system 
state-versus-time  graphs  can  be  updated  with  fixed  length  horizontal  line  increments. 
Although  the  clock  is  incremented  in  fixed  increments,  and  all  events  are  plotted  as  if 
they  occurred  at  the  end  of  At,  the  program  does  preserve  the  correct  chronological 
order  of  the  events  when  plotting. 

For  purposes  of  the  estimated  performance  measures,  the  program  keeps  track  of 
continuous  time,  and  does  not  round  event  occurrences  up  to  At.  In  this  sense,  the 
program  behaves  as  if  it  is  using  a  next-event  clock. 

B.  THE  COMPUTER 

The  computer  used  for  this  simulation  is  the  IBM  Personal  Computer  (IBM  PC), 
equipped  with  a  color  monitor.  In  order  to  run  the  simulation  the  IBM  PC  must  be 
configured  with  at  least  128  kilo-bytes  of  random  access  memory  (RAM),  one  5  1/4" 
floppy  disk  drive,  and  a  Color/Graphics  adapter.  The  color  monitor  and 
Color/Graphics  adapter  are  essential  since  the  program  uses  color  graphics. 

The  simulation  program,  JACKQUE  was  written  using  version  2.1  of  the 
PC-DOS  disk  operating  system.  Any  version  of  PC-DOS  earlier  than  2.1  will  also 
work.  Because  JACKQUE  directly  interacts  with  the  disk  operating  system,  there  can 
be  no  guarantee  that  the  simulation  will  run  on    IBM  Compatible  machines. 

C.  THE  PROGRAMMING  LANGUAGE 

JACKQUE  is  written  in  IBM  Advanced  Basic  (BASICA).  This  language  is 
provided  with  PC-DOS  and  possesses  all  of  the  graphics  capabilities  necessary  for  the 
simulation.    This  version  of  Basic  is  an  interpreter  language  and  runs  slower  than  a 
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compiled  language.     However  since  the  illusion  of  elapsed  time  is  desired  in  this 
particular  application,  slower  program  run  time  is  not  a  disadvantage. 

JACKQUE  also  uses  one  sub-routine,  called  RNGEN,  which  is  written  in  8088 
assembly  language. 

D.  PROGRAM  STRUCTURE 

The  program  listing  for  JACKQUE  is  contained  in  Appendix  A.  The  program,  as 
it  appears  in  the  Appendix,  has  been  extensively  edited  with  comments  to  assist  the 
interested  reader  in  understanding  the  program.  The  program  contained  on  the 
distribution  diskette  does  not  include  the  comments,  since  the  inclusion  of  the 
comments  exceeds  the  memory  available  within  the  Basic  stack.  Additionally,  extensive 
use  of  subroutines  were  used  in  an  effort  to  facilitate  following  the  program  flow.  The 
program  consists  of  five  major  modules.  These  modules  are  labeled  Main  Program 
Module  (MPM),  Clock  Module  (CM),  Arrival  Module  (AM),  Service  Center  1  Module 
(SCIM),  and  Service  Center  2  Module  (SC2M).  Any  subroutine  called  exclusively  by 
one  of  the  five  modules  is  physically  located  within  that  module.  After  the  fifth 
module,  all  subroutines  that  are  called  from  more  than  one  module  are  listed,  and 
labeled  accordingly.  All  subroutines  have  been  labeled  with  appropriate  titles.  After 
each  subroutine  title  there  is  a  slash  (/),  followed  by  the  abbreviation  of  the  module  or 
modules  from  which  it  is  called. 

E.  DESCRIPTION  OF  GRAPHICS  DISPLAY 

There  are  two  configurations  of  graphics  available  to  the  user.  When  the 
simulation  is  first  started,  the  graphics  display  screen  will  appear  as  seen  in  Figure  3.1. 
The  upper  half  of  the  screen  contains  a  pictorial  representation  of  the  status  of  the 
Jackson  network  for  a  particular  point  in  time.  A  circle  in  a  box  represents  a  person 
being  served  at  that  service  center.  A  circle  in  front  of  a  box  represents  a  person 
waiting  at  that  service  center.  There  may  be  more  customers  waiting  than  can  be 
physically  displayed  on  the  screen.  The  number  contained  in  the  circle  indicates  the 
number  of  times  that  the  particular  customer  has  looped  back  to  service  center  1  from 
service  center  2.  The  absence  of  a  number  indicates  that  the  customer  is  just  entering 
the  network  and  has  not  yet  experienced  any  feedback. 


"There  are  compiled  versions  of  Basic  in  existence 
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Figure  3.1     Screen  Display  with  Statistics  Table. 

In  Figure  3.1,  there  are  twelve  customers  waiting  at  service  center  1  as  indicated 
by  Qj  =  12,  in  addition  to  the  one  customer  currently  receiving  service.  At  service 
center  2  there  are  four  customers  waitmg  for  service  in  addition  to  the  one  customer 
receiving  service.  Of  the  customers  visible,  four  have  experienced  feedback  twice  while 
one  has  experienced  feedback  once. 

The  bottom  left  corner  of  the  screen  contains  the  state-versus-time  graph  for 
service  center  1.  The  system  state  values  for  the  service  center  are  displayed  on  the 
ordinate  of  the  graph,  and  32  increments  representing  time  on  the  abscissa.  The 
display  shows  that  at  time  =  3818,  the  state  of  service  center  1  jumped  from 
state  8  to  7. 

The  table  in  the  lower  right  corner  contains  the  estimated  and  theoretical 
quantities  discussed  in  Section  II-D.  The  variable  names  appearing  in  this  table  are 
consistent  with  those  defined  in  Chapter  II. 

Once  the  simulation  is  running,  the  user  can  change  the  configuration  of  the 
graphics  screen  to  that  shown  in  Figure  3.2.  In  this  case,  the  state-versus-time  graph 
for  service  center  2  is  displayed  in  the  lower  right  corner  of  the  screen.  If  the 
state-versus-time  graphs  for  both  service  centers  were  superimposed  on  each  other,  the 
resulting  graph  would  be  the  overall  system  state-versus-time  graph. 
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Figure  3.2     Screen  Display  with  out  Statistics  Table. 

A  more  complete  description  of  the  graphic  displays,  and  user  commands 
afTecting  them,  are  contained  in  Chapter  IV,  the  user's  manual. 

F.       PROGRAM  DESCRIPTION 

The  main  program  module  contains  numerous  subroutines  that  accomplish  most 

of  the  routine  housekeeping  chores.   These  include: 

Checking  to  see  if  BASICA  has  been  loaded  from  DOS,  and  if  the  machine  is 
configured  with  a  color/ graphics  adapter. 

Printing  the  title  screen. 

Loading  user  defined  functions. 

Loading  the  random  number  generator  sub-routine. 

Initializing  and  dimensioning  variables. 

Printing  the  main  program  menu. 

Creating  symbols  that  will  be  used  latter  in  the  graphic  display. 

Generating  the  graphics  display. 

Computing  the  theoretical  values  of  the  tabulated  quantities. 

The  Program  Menu  presents  the  user  with  a  number  of  options:  (1)  view  the 
program  instructions,  (2)  select  a  new  seed  for  the  random  number  generator,  (3)  start 
the  simulation  using  the  default  parameters,  (4)  change  the  model  parameters,  and  (5) 
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end  the  program.  The  on-screen  instructions  provide  a  condensed  version  of  the  user 
instructions  found  in  Chapter  IV,  the  user's  manual.  The  main  program  module  also 
contains  subroutines  to  perform  the  tasks  selected  from  the  program  menu.  The  last 
task,  performed  by  the  main  program  module  is  to  transfer  control  of  the  program  to 
the  clock  module. 

The  programming  for  some  of  the  housekeeping  and  administrative  tasks  found 
in  the  main  program  module  are  taken  in  whole  or  in  part  from  the  program  written  by 
R.E.  Nelsen  [Ref  2].  In  particular  the  subroutine  to  check  for  BASICA  and  the 
graphics  card,  the  subroutine  to  load  the  assembly  language  random  number  generator 
program,  and  the  subroutine  to  change  the  seed  of  the  random  number  generator  were 
taken  essentially  without  change. 

Once  a  simulation  run  is  initialized  and  running,  the  clock  module  controls  the 
How  of  the  program  for  the  remainder  of  the  simulation  run.  The  first  function  of  the 
clock  module  is  to  generate  an  initial  exogenous  customer  arrival  time  and  determine  if 
that  customer  will  exit  after  completing  service  at  service  centers  1  and  2.  This  is  a  one 
time  initialization. 

The  clock  module  determines  the  next-event  time  as  being  the  smallest  of  the 
next  arrival  time,  next  service  completion  time  at  service  center  2,  or  next  service 
completion  time  at  service  center  2.  The  next-event  time  thus  determined  is  always 
rounded  up  to  an  integer  value  so  to  be  compatible  with  the  fixed  increment  simulation 
clock. 

Every  time  the  simulation  clock  advances  one  time  increment,  the  clock  module 
updates  the  state-versus-time  graphs  for  those  service  centers  selected  to  be  displayed 
by  the  user.  Recall  that  the  state  of  a  particular  service  center  is  the  total  number 
waiting  and  being  served  by  that  service  center  at  a  particular  point  in  time.  When  the 
elapsed  time  reaches  32  time  units,  the  plotted  lines  on  these  graphs  will  be  at  the  right 
most  edge  of  the  graphs.  When  this  condition  occurs,  the  clock  module  shifts  the 
displayed  state  graphs  one  time  unit  to  the  left,  and  then  plots  the  new  time  increment. 
This  shifting  of  the  graphs  by  one  time  unit  to  the  left  continues  for  all  successive  time 
units.  In  this  manner  the  service  center  state-versus-time  graphs  always  contain  the 
most  recent  32  time  units  of  information,  while  more  dated  information  is  permanently 
lost.  When  the  simulation  clock  eventually  equals  the  next-event  time,  the  clock 
module  shifts  control  to  either  the  arrival  module,  service  center  one  module,  or  service 
center  2  module  as  appropriate. 
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Every  time  the  clock  module  increments  the  simulation  clock,  it  also  checks  the 
keyboard  for  certain  commands  issued  by  the  user.  By  striking  the  appropriate 
keyboard  keys,  the  user  can  suppress  the  audible  tones  associated  with  movement  in 
the  network,  pause  the  simulation,  continue  the  simulation,  change  the  display 
configuration,  speed  up  the  simulation,  and  end  the  simulation.  These  user  commands 
are  described  in  detail  in  Chapter  IV,  the  user's  manual. 

There  are  only  five  possible  events  that  can  happen  in  the  simulation: 

1.  An  exogenous  customer  can  enter  the  network. 

2.  A  customer  can  complete  service  at  service  center  1  and  exit  the  network. 

3.  A    customer    can    complete    service    at    service    center    1    and    proceed    to 

service  center  2. 

4.  A  customer  can  complete  service  at  service  center  2  and  exit  the  network. 

5.  A  customer  can  complete  service  at  service  center  2  and  loop  back  to  service 
center  1. 

Event  1  is  handled  by  the  arrival  module.    Events  2  and  3  are  handled  by  the  service 

center  1  module.   Events  4  and  5  are  handled  by  the  service  center  2  module. 

The  arrival  module  is  called  when  an  exogenous  customer  enters  the  network.  It 
flashes  the  arrow  leading  into  service  center  1.  If  no  one  is  presently  in  service,  a  blue 
ball  representing  the  customer  is  placed  in  the  square  box  corresponding  to  service 
center  1.  If  there  is  a  customer  already  in  service  but  the  waiting  queue  has  less  than 
eight  people  in  it,  a  blue  ball  is  placed  at  the  end  of  the  queue  and  the  queue  counter, 
Qj,  is  incremented  by  one.   Otherwise  only  Q^  is  incremented. 

If  the  exogenous  customer  immediately  receives  service,  rather  than  entering  the 
queue,  the  arrival  module  calls  the  random  number  generator  and  computes  the  length 
of  service.  In  all  cases  the  arrival  module  calls  the  random  number  generator  and 
computes  the  inter-arrival  time  of  the  next  exogenous  customer.  The  inter-arrival  time, 
and  service  time  if  applicable,  are  added  to  the  current  time  (continuous  not  rounded 
time),  to  determine  when  the  service  will  be  completed  and  when  the  next  exogenous 
arrival  will  take  place.  When  the  program  returns  to  the  clock  module,  these  newly 
created  event  times  will  be  compared  to  the  service  completion  time  at  service  center  2 
to  determine  what  event  should  occur  next  in  the  simulation. 

The  service  center  1  module  is  called  when  the  next  scheduled  event  is  a  service 
completion  at  service  center  I.  The  departing  customer  can  either  depart  the  network 
or  continue  on  to  service  center  2.  In  either  case,  the  module  erases  the  departing  ball, 
shifts  any  customers  waiting  for  service  center  1  to  the  right,  and  decrements  Q^^  by  one 
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as  appropriate.  If  there  is  now  a  customer  receiving  service  at  service  center  1,  the 
random  number  generator  is  called  and  a  new  service  completion  time  is  generated. 

If  the  customer  departs  the  system,  the  appropriate  arrow  is  flashed.  If  the 
customer  proceeds  to  service  center  2  the  arrow  leading  from  service  center  1  to  service 
center  2  is  flashed.  If  no  one  is  presently  in  service  at  service  center  2,  the  transferred 
ball  will  be  placed  in  the  square  box  representing  service  center  2.  If  there  is  a 
customer  already  in  service,  but  the  waiting  queue  in  front  of  service  center  2  has  less 
than  five  people  in  it,  the  ball  is  placed  at  the  end  of  the  queue,  and  the  queue  counter 
Q  ,  is  incremented  by  one.   Otherwise,  only  Q^  will  be  incremented  by  one. 

If  the  transferred  customer  immediately  receives  service  at  service  center  2,  rather 
than  exiting  or  entering  the  queue  in  front  of  service  center  2,  the  service  center  2 
module  calls  the  random  number  generator  and  determines  a  service  completion  time 
for  that  customer.  The  service  completion  time  is  added  to  the  current  time  to 
determine  when  the  service  will  be  completed.  When  the  program  returns  to  the  clock 
module,  this  newly  created  event-time  will  be  compared  to  the  next  exogenous  arrival 
time  and  service  completion  time  at  service  center  2,  to  determine  what  event  should 
occur  next  in  the  simulation.  Before  returning  to  the  clock  module,  the  service  center  1 
module  calls  the  random  number  generator  and  determines  whether  or  not  the  next 
customer  leaving  service  center  1  exits  the  network  or  proceeds  on  to  service  center  2. 

The  service  center  2  module  is  called  when  the  next  scheduled  event,  as 
determined  by  the  clock  module,  is  a  service  completion  at  service  center  2.  The 
departing  customer  can  either  exit  the  network  or  loop  back  to  service  center  1.  In 
either  case  the  service  center  2  module  erases  the  departing  ball,  shifts  any  customers 
waiting  for  service  center  2  to  the  right,  and  decrements  Q^  as  appropriate.  If  there  is 
now  a  customer  receiving  service  at  service  center  2,  the  random  number  generator  is 
called  and  a  new  service  completion  time  is  generated. 

If  the  customer  departs  the  system,  the  service  center  2  module  flashes  the 
appropriate  arrow.  If  the  customer  loops  back  to  service  center  1,  the  bottom  arrow 
from  service  center  two  to  service  center  one  is  flashed.  If  no  customer  is  presently  in 
service  at  service  center  1,  the  transferred  ball  is  placed  in  the  square  box  representing 
service  center  1.  If  there  is  a  customer  already  in  service,  but  the  waiting  queue  has 
less  than  8  people  in  it,  the  transferred  ball  is  placed  at  the  end  of  the  queue,  and  the 
queue  counter,  Q^^,  is  incremented  by  one.   Otherwise,  only  Qj^  is  incremented. 
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If  the  transferee!  customer  immediately  receives  service  at  service  center  1,  rather 
than  exiting  or  entering  the  queue  in  front  of  service  center  1,  the  service  center  2 
module  calls  the  random  number  generator  and  determines  a  service  completion  time 
for  that  customer.  The  service  completion  time  is  added  to  the  current  time  to 
determine  when  the  service  will  be  completed.  When  the  program  returns  to  the  clock 
module,  this  newly  created  event  time  will  be  compared  to  the  next  exogenous  arrival 
time  and  service  completion  time  at  service  center  1,  to  determine  what  event  should 
occur  next.  Before  returning  to  the  clock  module,  the  service  center  2  module  calls  the 
random  number  generator  and  determines  whether  or  not  the  next  customer  leaving 
service  center  2  exits  the  network  or  loops  back  to  service  center  1. 

The  service  center  2  module  performs  one  other  important  function.  It  keeps 
track  of  how  many  times  a  customer  loops  back  to  service  center  1,  allowing  this 
number  to  be  physically  placed  inside  the  balls  representing  customers.  When  a 
customer  first  arrives  and  has  had  no  opportunity  to  loop  from  service  center  2  back  to 
service  center  1,  he  is  depicted  as  a  solid  blue  ball  with  no  number  printed  on  it.  If  a 
customer  has  looped  from  service  center  2  to  service  center  1  exactly  once,  he  will  be 
depicted  as  a  purple  ball  with  a  number  "1"  printed  in  the  center  of  it.  If  a  customer 
has  looped  exactly  twice,  he  will  be  depicted  as  a  purple  ball  with  a  number  "1"  printed 
in  the  center  of  it,  etc. 

All  of  the  event  modules  create  audible  sounds  that  correspond  to  customer 
motion  in  the  network.  There  are  only  two  distinct  tones.  A  short  high  pitched  tone 
to  indicate  a  departure  from  the  network,  and  a  longer  low  pitched  tone  to  indicate 
arrivals  from  both  outside  the  network  and  arrivals  from  another  service  center. 

When  a  ball  arrives  at  a  service  center  and  finds  the  visible  queue  space  at 
capacity  (  8  for  service  center  1  and  5  for  service  center  2),  the  modules  store  these 
customers  in  a  queue  that  is  not  visible  to  the  user.  These  invisible  queues,  one  for  each 
service  center,  preserve  the  order  and  identity  of  the  overflow  customers.  They  are 
overflow  only  in  the  sense  that  there  is  not  enough  room  on  the  computer  screen  to 
graphically  display  them  all.  If  the  queue  counters  (Q^^,  Q^)  do  not  agree  with  what  is 
physically  shown  on  the  screen,  it  is  because  there  are  customers  in  the  invisible  queue 
or  queues.   The  queue  counters  are  correct. 

After  returning  to  the  clock  module,  the  simulation  updates  the  estimates  of  the 
tabulated  performance  values,  using  information  generated  in  the  event  modules,  and 
prints  them  to  the  screen.    As  previously  noted,  the  clock  module  also  updates  the 
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service  center  state-versus-time  graphs.  Each  of  the  event  modules  records  the 
information  necessary  for  the  clock  module  to  know  where  vertically  on  the  graphs  it 
should  draw  the  horizontal  line  segments  representing  one  time  unit.  These  graphs, 
because  of  space  constraints,  can  not  plot  state  variables  greater  than  nine.  If  the  state 
variable  is  equal  to  9,  the  clock  module  plots  a  purple  tick  mark  at  a  vertical  position 
corresponding  to  9.  For  state  variables  greater  than  9,  the  clock  module  plots  a  blue 
tick  mark  also  at  a  vertical  position  corresponding  to  9.  Therefore,  observing  a  blue 
tick  mark  only  indicates  that  the  service  center  state  is  10  or  greater. 

G.       PROGRAiMMING  CONSIDERATIONS 

At  all  times  during  the  simulation,  JACKQUE  maintains  a  next-event  time  for 
each  of  the  three  possible  events  (exogenous  arrivals,  and  service  completion  at  service 
centers  1  or  2).  The  clock  module  than  takes  the  minimum  of  these  three  times  to 
determine  the  next  event.  Because  of  the  structure  of  the  network,  there  are  times 
when  an  event  should  never  be  allowed  to  occur. 

For  example,  if  no  customer  is  currently  being  served  at  a  service  center,  the 
program  must  be  prevented  from  selecting  the  service  completion  time  of  that  service 
center  as  the  minimum  of  the  next-event  times.  If  the  program  did  select  the  service 
completion  time  as  the  next-event,  the  simulation  would  attempt  to  move  a  customer 
from  the  service  center  when  in  fact  there  is  no  customer  there. 

To  prevent  the  selection  of  an  impossible  event,  the  program  sets  the  next-event 
time  of  all  impossible  events  equal  to  10^^,  which  for  all  practical  purposes  is  equal  to 
infinity.  The  clock  module  takes  the  minimum  of  the  next-event  times.  It  will  never 
select  10^^,  because  the  next  exogenous  arrival  time  will  always  be  less  than  this.  In 
this  particular  model,  one  or  both  of  the  next  service  times  may  periodically  be  set  to 
infinity. 

When  the  clock  module  selects  the  minimum  of  the  three  next-event  times,  there 
is  a  remote  possibility  that  there  will  be  a  two-way,  or  even  three-way  tie.  In  the 
unlikely  event  that  a  tie  does  occur,  the  program  decides  randomly  which  of  the  events 
should  occur  first. 

Because  of  the  discrete  nature  of  the  graphics  display,  events  can  only  be 
depicted  as  having  occurred  at  the  end  of  a  discrete  time  tick.  For  this  reason,  within 
one  time  nek,  At,  many  events  may  actually  occur.  Some  of  these  events  may  occur  at 
the  very  begining  of  At,  some  in  in  the  middle,  and  perhaps  some  near  the  end  of  At. 
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Although  their  chronological  order  of  occurrence  is  preserved,  they  are  all  plotted 
sequentially  as  if  they  occurred  at  the  end  of  At.  For  this  reason  what  is  observed  on 
the  screen  is  not  an  absolutely  accurate  representation  of  a  Jackson  network  with 
exponentially  distributed  inter-arrival  times  and  service  times.  What  is  graphically 
displayed  is  a  Jackson  network  with  geometrically  distributed  inter-arrival  and  service 
times,  the  geometric  distribution  being  the  discrete  analog  of  the  continuous 
exponential  distribution.  Despite  this  fact,  the  simulation  still  produces  a  realistic 
visual  impression  of  a  true  Jackson  network. 

For  a  discrete  event  stochastic  process,  the  state  of  the  process  can  only  change 
by  discrete  jumps  of  plus  or  minus  one.  However,  because  JACKQUE  considers  all 
events  occurring  in  At  to  have  occurred  at  the  end  of  At,  the  state-versus-time  graphs 
might  exhibit  jumps  greater  than  plus  or  minus  one.  For  example  if  service  center  1 
experiences  3  arrivals  and  1  service  completion,  all  during  At,  the  state-versus  time 
graph  for  service  center  1  will  jump  2  units  when  the  horizontal  line  segment 
corresponding  to  At  is  plotted  (see  Figure  3.3).  Similarly,  the  state-versus-time  graph 
might  plot  the  next  tick  mark  at  the  same  vertical  height  as  the  previously  plotted  tick 
mark,  indicating  no  change  when  in  fact  there  had  been  a  change  during  At.  For 
example,  if  service  center  1  experiences  one  arrival  and  one  service  completion  during 
At,  there  is  no  net  change  and  the  time  tick  will  be  plotted  at  the  same  vertical  height 
as  before  (see  Figure  3.4). 

This  undesireable  effect  can  be  minimized  by  reducing  the  probability  of  having 
multiple  events  occur  during  At.  This  can  be  accomplished  by  either  reducing  the  size 
of  At  or  through  the  judicious  selection  of  the  network  parameters,  as  discussed  in 
Section  IV-I. 

H.       GENERATION  AND  USE  OF  RANDOM  NUMBERS. 

As  demonstrated  by  Nelsen  [Ref.  2],  the  random  number  generator  residing  in 
BASICA  produces  pseudo-uniform  (0,1)  numbers  with  significant  lag-1  correlation. 
This  makes  it  a  poor  random  number  generator  and  unsuitable  for  use  in  simulation. 

JACKQUE  uses  a  random  number  generator  modified  by  Nelsen,  and  originally 
programmed  by  B.O.  Shubert  of  the  Naval  Postgraduate  School  for  use  on  the 
IBM-PC.  This  random  number  generator,  named  RNGEN,  is  written  in 
assembly  language. 
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RNGEN  is  a  Prime  Multiplicative  Linear  Congruential  Generator  (PMMLCG) 
based  on  ihe  following  recursive  formula: 

X-  +  ^=  75xMod(2^^-l),  (eqn3.1) 

where  X-  is  the  old  seed  and  X-^  j^  is  the  new  seed.  A  uniform  (0,1)  random  number  is 
obtained  by  dividing  the  new  seed  by  the  modulus.  This  particular  PMMLCG  is 
considered  to  exhibit  good  statistical  behavior  and  is  currently  used  in  the 
LLRANDOM  package,  IBM's  GGL,  and  IMSL's  GGLBS  [Ref  6:  p.226]. 

JACKQUE  uses  Uniform  (0,1)  random  variables  for  three  different  purposes: 

1.  To  generate  exponentially  distributed  random  variables. 

2.  To    determine   whether   a    customer   exits    or   remains   in   the    network   after 
completing  service  at  a  service  center. 

3.  To  break  ties  in  the  event  two  or  more  event  times  are  minimum. 

The  inter-arrival  times  and  service  times  of  customers  in  a  Jackson  network  are 
assumed  to  be  exponentially  distributed.  To  obtain  these  exponential  random 
variables.  Uniform  (0,1)  random  variables  produced  by  RNGEN  are  transformed 
through  the  inverse  probabiUty  method.  Let  U  be  a  uniformly  distributed  (0,1)  random 
variable,  and  let  X  be  an  exponentially  distributed  random  number  with  parameter  X. 
The  distribution  function  for  X  is: 

F(x)  =   1  -  e-^'^   ,x  >  0.  (eqn  3.2) 

An  exponential  random  number  is  generated  by  setting  u  =  F(x)  and  solving  for  x: 

u  =    1  -  e"^'^  (eqn  3.3) 


1-u  =  e'^^  (eqn  3.4) 


In(l-u)  =  -Xx  (eqn  3.5) 

x  =  -l/}iln(l-u)  (eqn  3.6) 
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X  =  -1/Xln(u).  (eqn  3.7) 

It  is  possible  to  replace  (1-u)  with  u  because  both  (1-u)  and  u  have  the  same  Uniform 
(0,1)  distribution. 

When  a  customer  completes  service  at  a  service  center,  he  can  either  exit  the 
system  or  go  to  the  next  service  center.  To  determine  which  occurs  JACKQUE 
compares  a  Uniform  (0,1)  random  number  to  the  fixed  probability  that  the  customer 
exits  after  completing  service  at  service  center  i  (P-).  If  U(0,1)  ^  Pj,  the  customer  is 
considered  to  have  exited.  Otherwise  the  customer  is  considered  to  have  proceeded  to 
the  next  service  center. 

A  similar  procedure  is  used  to  break,  a  tie  in  the  event  two  or  three  event-times 
are  simultaneously  minimum.  In  the  case  of  a  two-way  tie,  event  A  is  chosen  if 
U(0,1)  <  1/2,  and  event  B  is  chosen  otherwise.  In  the  case  of  a  three-way  tie,  event  A 
is  chosen  if  0  <  U(0,I)  <  1/3,  event  B  is  chosen  if  1/3  <  U(0,1)  <  2/3,  and  event  C  is 
chosen  otherwise.  Because  all  events  do  not  have  equal  probabilities  of  occurring, 
using  equal  probabilities  is  not  totally  accurate,  but  considered  sufficient  for  the 
purposes  of  this  simulation. 
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A.  INTRODUCTION 

JACKQUE  is  a  discrete  event  stochastic  simulation  of  a  Jackson  network.  A 
Jackson  network  is  a  network  of  queues  where  a  customer  may  receive  service  from 
numerous  service  centers,  possibly  visiting  the  same  service  center  more  than  once, 
before  exiting  the  network.  The  program  produces  a  variety  of  visual  displays  on  the 
computer  screen  intended  to  provide  the  user  with  a  better  understanding  of  the 
dynamic  nature  of  the  model. 

This  manual  is  intended  to  assist  the  user  in  the  actual  use  and  operation  of  the 
program.  In  describing  the  user  interaction  with  the  program,  the  following  notation  is 
used: 

•  When  a  command  is  to  be  entered  (typed)  exactly  as  shown,  it  will  be  enclosed  in 
single  quotes  ('  '). 

•  When  a  specific  key  on  the  keyboard  is  to  be  pressed,  the  letter  or  letters  on  the 
key  will  be  encased  by  brackets  (  <   >  ). 

For  example,  'Load  "JACKQUE. BAS",r'  <Enter>  ,  means  type  the  command  exactly 

as  it  appears  between  the  single  quotes  and  then  press  the  ENTER  key.  PC-DOS  and 

JACKQUE  make  no  distinction  between  upper  and  lower  case  letters,  so  either  may  be 

used. 

B.  MODEL  DESCRIPTION 

The  Jackson  network  modeled  by  this  simulation  consists  of  two  service  centers, 
each  containing  one  server.  An  arriving  customer  from  outside  the  network  must  first 
go  to  service  center  1  where  he  enters  into  a  First  in  First  Out  (FIFO)  queue.  Upon 
completion  of  service  at  service  center  1,  the  customer  either  exits  the  network  with 
probability  P^,  or  proceeds  on  to  service  center  2  with  probability  l-P^^,  where  he  enters 
another  FIFO  queue.  Upon  completion  of  service  at  service  center  2,  the  customer 
either  exits  the  network  with  probability  P  ,  or  loops  back  to  service  center  2  with 
probability  l-P-.  again  entering  a  FIFO  queue.  In  this  model,  a  customer  can  loop 
from  service  center  2  to  service  center  1  a  large  but  finite  number  of  times. 

There  are  several  assumptions  required  for  this  model.  Customers  from  outside 
the  network  are  assumed  to  arrive  at  service  center  1  according  to  a  Poisson  process 
with  rate  "k.  The  service  times  at  service  center  1  and  service  center  2  are  assumed 
exponentially  distributed  with  rates  yi^  and  \\.^  respectively.  When  a  customer  departs 
the  network  with  probability  P,  or  goes  to  the  next  service  center  with  probability  I-P, 
he  does  so  instantaneously.   These  probabilities  are  assumed  independent  of  the  history 
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of  the  network.  It  is  assumed  that  there  is  unlimited  waiting  space  at  both  service 
centers,  and  that  if  one  server  is  idle  he  does  not  help  a  busy  server  at  the  other  service 
center. 

For  this  model,  the  state  of  the  system,  N'(t),  is  the  sum  of  N^(t)  and  N^Ct).  Nj(t) 
(i=  1,2)  is  defined  as  the  number  of  customers  waiting  and  being  served  at  service 
center  i  at  time  t.   Therefore  N(t)  =  N^(t)  +  N^Ct)  =  n,  where  n  =  0,1,2,.... 

C.       HARDWARE  AND  SOFTWARE  REQUIREMENTS 

1 .  Hardware 

An  IBM  Personal  Computer  (IBM-PC)  with  at  least  128  Kilo-bytes  of  random 
access  memory  (RAM),  a  color/graphics  adapter  card,  a  color  monitor,  and  at  least 
one  5  1/4  inch  floppy  disk  drive  are  required  to  run  the  simulation.  The  color  monitor 
and  color/graphics  adapter  are  essential  since  the  program,  JACKQUE  ,  uses  color 
graphics.  The  128  K  of  memory  is  required  because  JACKQUE,  which  is  written  in 
BASICA,  uses  an  assembly  language  sub-routine  which  is  loaded  outside  of  Basic's 
64  Kilo-byte  workspace. 

2.  Software 

The  primary  program  JACKQUE. BAS,  and  the  assembly  language 
sub-routine  RNGEN.SRT  are  provided  to  the  user  on  the  distribution  diskette.  In 
addition  to  these  programs,  the  user  must  provide  the  Disk  Operating  System  PC-DOS 
(2.1)  and  Advanced  Basic  (BASICA)  in  order  to  run  the  simulation.  Because 
JACKQUE  directly  interacts  with  the  Disk  Operating  System,  there  can  be  no 
guarantee  that  the  simulation  will  run  on  IBM  compatible  machines. 

3.  Program  Files 

The  following  program  files  are  included  on  the  distribution  diskette: 

•  JACKQUE. BAS  -  the  main  program  written  in  BASICA. 

•  RNGEN.SRT  -  a  random  number  generator  written  in  8088  assembly  language, 
used  as  a  sub-routine  to  the  main  program. 

•  JACKQUE.BAT  -  a  batch  file  that  loads  and  starts  the  simulation. 

•  AUTOEXEC.BAT  -  an  autoexec  batch  file  that  loads  and  starts  the  simulation 
immediately  after  the  system  is  started  or  restarted. 

The  first  two  files  are  absolutely  necessary  to  run  the  simulation.    The  two  batch  files 

are  two  different  methods  of  automatically  loading  and  starting  the  simulation.    The 

batch  files  are  not  essential,  but  enhance  user-friendliness.    A  source  code  listing  of 

RNGEN.SRT  is  contained  in  Reference  2. 
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D.      GETTING  STARTED 

1.  Making  an  Application  Diskette 

The  application  diskette  is  created  using  the  distribution  diskette.  After  being 
used  for  this  purpose,  the  distribution  diskette  should  be  stored  in  a  safe  place,  and  only 
used  to  create  additional  application  diskettes  as  needed. 

Although  there  are  many  different  ways  that  a  user  can  create  an  application 
diskette,  this  manual  will  describe  just  one.  A  self-starting  (bootable)  diskette 
containing  certain  operating  files,  the  BASICA  command  file  and  the  simulation 
program  files  is  considered  the  most  efiicient  format.  This  eliminates  the  requirement 
for  a  separate  diskette  containing  PC-DOS,  and  the  diskette  can  be  run  using  only  one 
disk  drive. 

To  make  a  bootable  diskette,  follow  the  instructions  in  the  PC-DOS  manual 
and  format  a  blank  diskette  using  the  FORMAT  command  with  the  /s  option.  The  /s 
option  copies  only  the  DOS  files  that  must  be  on  the  newly  formatted  disk.  These  are 
the  COMMAND.COM  file,  plus  several  hidden  files  that  are  necessary  for  DOS  to 
operate.  The  Advanced  Basic  file,  called  BASICA.COM,  also  located  on  the  DOS 
diskette  can  be  transferred  onto  the  application  diskette  using  the  copy  command.  The 
copy  command  is  also  used  to  transfer  the  simulation's  four  files  from  the  distribution 
diskette  onto  the  bootable  application  diskette.  The  user  now  has  a  totally  self 
contained,  self  starting  diskette  capable  of  running  the  simulation. 

2.  Starting  the  Simulation 

If  the  PC  is  turned  ofi',  place  the  bootable  diskette  into  drive  A  and  turn  on 
the  computer.-^  After  the  computer  does  it's  internal  checks  taking  nearly  60  seconds, 
the  computer  will  immediately  execute  the  DOS  commands  contained  in  the 
AUTOEXEC  batch  file.  These  DOS  commands  automatically  load  BASICA  and 
JACKQUE.BAS,  and  run  the  simulation. 

If  the  PC  is  already  turned  on  and  in  the  command  level  of  DOS,  place  the 
bootable  diskette  in  drive  A  and  type  'JACKQUE' <  ENTER>  .  The  computer  will 
immediately  execute  the  DOS  commands  contained  in  the  JACKQUE.BAT  file  and 
start  the  simulation.  Alternately,  with  the  bootable  diskette  in  drive  A,  the  computer 
can  be  restarted  by  simultaneously  striking  <Ctrl>  <alt>  <del>.  In  this  case  the 
simulation  is  started  using  the  AUTOEXEC.BAT  file. 


^This  section  assumes  that  drive  A  is  the  default  drive. 
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If  for  some  reason  the  computer  is  in  the  command  level  of  DOS  and  it  is  not 
desired  to  use  the  batch  files,  the  simulation  can  be  started  by  using  the  following 
sequence  of  commands: 

•  BASICA'  <RETURN> 

•  'Load  "d:JACKQUE.BAS",r'  <  ENTER>  where  d  is  the  label  of  the  drive  where 
JACKQUE.BAS  can  be  found. 

If  the  simulation  has  been  terminated  from  the  program  menu,  and  the 
computer  is  still  in  the  BASICA  environment,  another  simulation  run  can  be  initiated 
by  simply  typing  'RUN'  <  ENTER  >  .  If  you  have  returned  to  the  command  level  of 
DOS,  use  the  methods  previously  described. 

During  the  execution  of  the  program,  JACKQUE  attempts  to  load  the 
random  number  generator  sub-routine  from  the  diskette  in  the  A  disk  drive.  If  the 
program  fails  to  find  the  necessary  file  in  drive  A,  it  asks  the  user  for  the  correct  drive 
label  where  it  can  be  found.  Simply  respond  by  entering  the  proper  drive  label  and 
<  ENTER >.  If  it  can  not  find  the  sub-routine  on  the  specified  drive,  or  if  the 
specified  drive  is  not  valid,  the  program  will  again  ask  the  user  for  the  correct  drive 
label.  This  situation  will  never  arise  if  the  bootable  application  diskette  is  used  in 
drive  A. 

3.    Title  Screen 

Upon  starting  the  simulation,  the  first  thing  the  user  will  observe  on  the 
screen  will  be  the  title  screen  as  shown  in  Figure  4.1.  From  this  point  on,  the 
on-screen  prompts  will  guide  the  user  through  the  running  of  the  simulation. 

E.       PROGRAM  MENU 

Upon  leaving  the  title  screen,  there  will  be  a  few  second  delay  while  JACKQUE 
initializes  variables  and  loads  the  random  number  generator  sub-routine.  The  next 
screen  to  be  displayed  will  be  the  program  menu  as  shown  in  Figure  4.2.  The  options 
on  the  menu  allow  the  user  to  read  the  instructions,  select  a  new  starting  seed  for  the 
random  number  generator,  run  the  simulation  using  the  default  parameters,  run  the 
simulation  using  parameters  selected  by  the  user,  and  end  the  program. 

When  a  simulation  is  ended,  JACKQUE  always  returns  to  the  program  menu. 
This  allows  the  user  to  conduct  numerous  simulations  runs,  perhaps  with  different 
seeds  and/or  parameters  values,  during  one  running  of  JACKQUE. 

Upon  departing  the  menu,  and  before  running  the  first  simulation,  JACKQUE 
quickly  plots  and  then  erases  the  screen  shown  in  Figure  4.3.    Before  being  erased, 
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Figure  4.1     The  Title  Screen. 


PROGRAM  MENU 


< I >ns true t ions 

<N>ew  random  number  Generator  Seed 

<D>efault  Model  Parameters 

<  C> han ^e  Mo de 1  Par ame  t e rs 

<E>nd  the  Proyram 


Enter  your  selection: 


Figure  4.2    The  Program  Menu. 

these  symbols  are  stored  in  memory  and  are  used  by  JACKQUE  in  generating  graphics 
for  all  future  simulation  runs. 
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Figure  4.3     Graphic  Symbols. 

1.  Displaying  the  User  Insiruciions 

To  display  the  user  instructions  on  the  screen,  press  <  I  >  .  The  instructions 
displayed  are  a  condensed  version  of  what  is  contained  in  this  guide. 

2.  Selling  the  Random  Number  Generator  Seed 

A  new  seed  can  be  selected  by  pressing  <  N  >  .  The  program  will  prompt  you 
to  enter  the  value  of  the  seed.  The  seed  is  entered  by  selecting  a  number  m  the  range 
from  1  to  2,147,483.646  and  pressing  <  ENTER  >  .  If  a  non-integer  number  is  entered, 
it  will  be  truncated  to  an  integer. 

If  no  seed  is  selected,  a  default  starting  seed  is  used.  If  the  simulation  is 
stopped,  and  then  a  new  simulation  started  without  leaving  JACKQUE  and  without 
setting  a  new  seed,  the  generator  will  continue  as  if  the  old  simulation  had  never 
stopped.  That  is,  the  first  random  number  generated  in  the  new  simulation  is  the  same 
random  number  that  would  have  been  generated  next  in  the  old  simulation. 

By  using  the  same  seed  with  the  same  model  parameters,  the  simulation  can 
be  run  numerous  times,  obtaining  the  exact  same  results  each  time.  This  may  be  useful 
for  replicating  a  particular  model  behavior. 
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3.    Default  Model  Parameters 

If  this  option  is  selected  by  pressing  <D>,  the  simulation  immediately  starts 
with  the  model  parameters  set  to  the  following  default  values: 
•    X  =  .10 


• 

Ml  =  .20 

• 

^2=  .13 

• 

P^  =  .30 

• 

P,=  .5 

If  a  new  seed  is  desired  it  must  be  entered  before  striking  <  D>  . 

4.  Changing  the  Model  Parameters 

When  <C>  is  pressed,  the  user  is  prompted  for  values  of  the  arrival  rate, 
service  rate  at  service  center  1,  service  rate  at  service  center  2,  probability  a  customer 
exits  after  departing  service  center  1,  and  probability  a  customer  exits  after  departing 
service  center  2.  The  rates  can  be  any  positive,  non-zero  real  number.  The 
probabilities  must  be  between  zero  and  one.  The  quantities  are  entered  by  typing  in 
the  appropriate  value  and  then  pressing  <  ENTER >  in  response  to  each  of  the  five 
questions.  If  <ENTER>  is  struck  without  typing  in  a  number,  the  parameter  value 
will  be  considered  to  be  zero. 

For  a  Jackson  network  to  be  stable,  the  effective  arrival  rate  at  each  service 
center,  taking  into  account  arrivals  from  both  outside  the  network  and  from  other 
service  centers,  must  be  less  than  the  corresponding  service  rates.  If  this  is  not  the  case 
the  network  will  become  unstable  in  that  more  customers  are  joining  the  queue  than 
departing  it. 

In  the  case  of  an  unstable  queue,  there  are  no  steady  state  probabilities  and 
the  network  performance  measures  have  no  meaning.  The  stability  of  the  network  is 
controlled  by  the  interactions  of  the  five  model  parameters.  Setting  P^^  and  P^ 
simultaneously  equal  to  0  will  always  result  in  an  unstable  network.  Section  III-I 
discusses  the  selection  of  model  parameters.  If  the  selection  of  parameters  results  in  an 
unstable  network,  the  screen  shown  in  Figure  4.4  will  be  displayed. 

The  simulation  starts  immediately  after  entering  the  last  requested  parameter 
value.    If  a  new  seed  is  desired,  it  must  be  entered  before  striking  <  C>  . 

5.  Ending  the  Program 

To  end  the  program,  press  <  E>  .  This  will  clear  the  screen  and  leave  the  PC 
in  the  BASIC  command  mode.    It  also  returns  the  screen  to  the  standard  monochrome 
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Vou^  parameter*  selection  if»esults 
in  Sn  Unstable  (EXPLOSIUE>  queue. 
The  gr*aphical  pr»esen  tat  ion  oi  tne 
networ»k  will  still  be  accux^ate; 
Howevex*  the  network  statistics 
will  be  INULA  ID,  and  will  NOT  be 
displayed.  See  the  user*'  s  wanual 
for*  an  explanation. 

<C>ontinue  with  Siwulation 
<R>etur*n  to  the  Prograw  MENU 

ENTER  your  Selection 


Figure  4.4     Example  of  Unstable  Parameter  Warning. 

mode  with  column  length  equal  to  80.    If  it  is  desired  to  return  to  the  DOS  command 
level  type  system'  <  ENTER >  . 

F.       THE  GIL^PHICS  DISPLAY 

JACKQUE  uses  various  graphical  displays  to  allow  the  user  to  observe  the  flow 
of  customers  through  the  network.  When  the  simulation  is  first  started,  the  graphics 
display  will  appear  as  shown  in  Figure  4.5.  The  graphics  display  is  divided  into  three 
regions:  a  pictorial  display  of  the  network,  the  state-versus-time  graph  for  service 
center  1,  and  a  table  of  network  performance  measures. 

At  any  time  after  the  start  of  the  simulation,  the  user  can  chose  to  change  the 
configuration  of  the  graphics  display  to  that  shown  in  Figure  4.6.  The  table  of 
performance  measures  has  been  replaced  by  the  state-versus-time  graph  for  service 
center  2.  The  user  can  freely  change  the  configuration  back  and  forth  between  the  two 
modes  as  described  in  Section  IV-G.  The  following  descriptions  are  directed  at  the 
display  mode  shown  in  Figure  4.5. 

1.    The  Pictorial  Network  Display 

The  display  in  the  upper  half  of  the  screen  is  a  pictorial  representation  of  the 
Jackson  network  at  a  particular  instant  in  time.  The  left  square  box  represents  service 
center  1.   The  right  square  box  represents  service  center  2.   The  vertical  and  horizontal 
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Figure  4.5     Screen  Display  with  One  State  Diagram. 
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Figure  4.6     Screen  Display  with  Two  State  Diagrams. 

lines  connecting  these  two  boxes  represent  the  possible  paths  that  a  customer  can 
follow  through  the  network.    An  arrow  directed  into  a  service  center  is  the  path  of  an 
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arriving  customer.  An  arrow  directed  out  of  a  service  center  whose  head  is 
unconnected,  is  the  path  of  a  customer  departing  the  network.  Arrows  connecting  the 
two  service  centers  are  the  paths  of  customers  moving  from  one  service  center  to  the 
other.  The  lowest  horizontal  line  directed  from  right  to  left  is  the  feedback  loop  from 
service  center  2  to  service  center  1. 

Customers  are  represented  as  colored  solid  balls.  A  customer  arriving  from 
outside  the  network,  for  the  first  time  is  represented  as  a  solid  blue  ball.  If  a  customer 
loops  from  service  center  2  to  service  center  1,  the  color  of  the  ball  is  changed  to 
purple,  and  a  number  "1"  is  placed  inside  the  ball.  Each  time  this  customer  loops  back 
from  service  center  2  to  service  center  i  the  number  inside  the  purple  ball  is 
incremented  by  one,  up  to  a  maximum  of  nine.  Thus  a  purple  ball  with  a  number  "5" 
in  it,  represents  a  customer  who  has  been  fed  back  from  service  center  2  to  service 
center  1  five  times.  When  a  customer  exits  the  network,  the  corresponding  ball 
disappears  from  the  screen. 

A  customer  actually  receiving  service  from  a  service  center  is  placed  inside  the 
appropriate  square  box.  If  a  customer  is  waiting  for  service  at  a  service  center,  he  is 
placed  in  the  waiting  Hne  immediately  to  the  left  of  the  appropriate  service  center 
(square  box).  The  counters  under  the  two  waiting  lines,  Q^  and  Q^,  display  the 
number  of  customers  waiting  for  the  respective  service  centers.  These  counters  do  not 
include  the  person  actually  receiving  service.  Because  of  the  Hmited  size  of  the 
computer  screen,  it  is  possible  for  the  values  of  Q^  and  Q^  to  be  greater  than  the 
number  of  customers  actually  shown  waiting  as  balls.  These  invisible  waiting  customers 
are  still  accounted  for,  and  will  enter  the  visible  waiting  line  when  space  permits. 

The  five  model  parameters,  X,  ]i^,  \i^,  P^^,  and  P^  are  displayed  in  appropriate 
places  on  the  pictorial  network  display.  The  current  time  of  the  system  is  printed  at 
the  bottom  of  this  display. 

Each  time  a  customer  (ball)  moves  in  the  network,  the  motion  is  represented 
by  moving  the  affected  balls,  flashing  the  appropriate  path  arrow,  and  sounding  a 
distinctive  tone.  There  are  only  two  distinctive  tones.  A  short  high  pitched  tone  to 
indicate  a  departure  from  the  network,  and  a  long  low  pitched  tone  to  indicate  arrivals. 
Arrivals  consists  of  both  arrivals  from  outside  the  network  and  arrivals  from  other 
service  centers. 

In  Figure  4.5,  both  service  centers  are  busy  with  four  customers  waiting  for 
service  center  2  {Qy=4),  and  twelve  customers  waiting  for  service  center  1  (Qj^=  12) 
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although  only  eight  are  physically  shown.  Of  the  visible  customers,  one  has 
experienced  feedback  twice,  four  have  experienced  feedback  once,  and  nine  are  passing 
through  the  network  for  the  first  time. 

2.  The  State- Versus-Time  Graph 

The  lower  left  corner  of  the  graphics  display  contains  the  state-versus-time 
graph  for  service  center  1.  The  values  of  the  state  variable  N(t),  the  number  of 
customers  waiting  and  being  served  at  service  center  1  at  time  t,  are  plotted  on  the 
ordinate.  Time  is  plotted  on  the  abscissa,  and  is  represented  as  32  discrete  divisions. 
Each  of  these  time  segments  corresponds  to  one  time  tick  of  the  simulation  clock. 

After  each  increment  of  the  simulation  clock,  the  program  determines  the 
current  state  of  service  center  I.  The  program  then  uses  this  information  to  plot  a 
horizontal  line  one  time  tick  in  length  at  the  appropriate  height  on  the 
state-versus-time  graph.  In  Figure  4.5,  at  time  3818,  the  state  of  service  center  1  went 
from  8  to  7. 

When  the  elapsed  time  reaches  32  time  units,  the  plotted  line  on  this  graph 
will  be  at  the  right  most  edge  of  the  graph.  When  this  condition  occurs,  the  program 
shifts  the  entire  graph  one  time  unit  to  the  left,  and  plots  the  new  time  increment  at  the 
appropriate  height.  In  this  manner  the  graph  always  contains  the  most  recent  32  time 
units  of  information  on  the  state  of  service  center  1.  An  identical  description  pertains 
to  the  state-versus-time  graph  for  service  center  2,  which  can  be  displayed  at  the  users 
option. 

The  state-versus-time  graphs,  because  of  space  constraints,  can  not  plot  state 
variables  greater  than  9.  If  the  state  variable  is  equal  to  9,  a  purple  tick  mark  at  a 
vertical  position  corresponding  to  9  is  plotted.  For  state  variables  with  values  greater 
than  nine,  a  blue  tick  mark  is  plotted,  also  at  a  vertical  position  corresponding  to  9. 
Therefore  observing  a  blue  tick  mark  only  indicates  that  the  service  center  state  is  10  or 
greater. 

3.  Network  Performance  Measures 

The  lower  right  hand  corner  of  the  graphics  display  contains  a  table  of 
estimated  and  theoretical  equiUbrium  quantities.  These  quantities  are  just  a  few  of  the 
many  that  can  be  used  to  evaluate  the  performance  characteristics  of  the  network.  The 
variable  names  used  in  this  table  are  defined  as: 

R-  =  the  effective  arrival  rate  of  customers  at  service  center  i  (i=  1,2). 
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S;  =  the  expected  number  of  times  a  customer  will  pass  through  service 
^       center  i  before  exiting  the  network  (i=  1,2). 

L:  =  the  expected  number  of  customers  both  waiting  and  being 
served  at  service  center  i  (i=  1,2). 

L  =  the  expected  number  of  customers  in  the  network,  both  waiting 
and  being  served. 

W:  =  the  expected  amount  of  time  a  customer  spends  waiting  and  being 
served  at  service  center  i  (i=  1,2). 

W  =  the  total  expected  amount  of  time  a  customer  spends  in  the  network,  both 
waiting  and  being  served. 

Z  =  the  proportion  of  time  service  center  1  is  busy  and  service 
center  z  is  idle. 

The  quantities  that  appear  under  the  column  header  Limit,  are  the  limiting 
values  of  these  quantities.  They  are  computed  using  results  from  Little's  Law  and  from 
the  analysis  of  simple  M/M/1  queues.  They  are  the  theoretical  values  that  would  be 
obtained  if  the  network  reached  steady  state. 

The  quantities  that  appear  under  the  column  header  Estimate  are  estimates  of 
these  theoretical  Umiting  values.  They  are  computed  with  data  collected  from  the 
simulation.  These  estimates  are  updated  after  every  event  occurrence,  where  events  are 
defined  to  be  the  arrival  of  a  customer  from  outside  the  network  and  service 
completions  at  either  service  center. 

At  certain  points  during  the  simulation,  some  of  the  estimates  might  be 
negative  and/or  extremely  small  or  large  in  magnitude.  If  an  estimate  exceeds  the 
format  capabiUties  of  JACKQUE  by  being  too  large  or  small,  the  estimate  will  be 
displayed  on  the  screen  as  "9999.99".  Some  of  the  estimates  will  temporarily  become 
negative  if  the  estimated  effective  arrival  rate  at  a  service  center  is  greater  than  the 
estimated  service  rate  of  the  service  center.  To  avoid  confusion,  the  absolute  value  of  a 
negative  estimate  will   be  displayed  on  the  screen. 

Because  of  the  manner  in  which  L^^,  L^,  W^^,  W^,  and  Z  are  estimated,  the 
estimates  might  show  considerable  variance  around  the  theoretical  value.  This 
variability  will  decrease  as  the  simulation  progresses  in  time. 

If  in  the  selection  of  model  parameters,  an  effective  arrival  rate  is  greater  than 
the  service  rate  at  the  applicable  service  center,  the  user  is  given  the  option  to 
terminate  the  current  simulation  or  to  continue  as  described  in  Section  IV-E-4.    If  the 
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user  opts  to  continue  with  the  present  model  parameters,  the  resulting  theoretical  and 
estimated  values  will  no  longer  be  valid.  In  this  case  the  graphics  display  will  appear 
as  seen  in  Fieure  4.7. 
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Figure  4.7    Example  of  Graphics  Screen  for  Unstable  Network. 


G.       KEYBOARD  COMMANDS 

While  the  simulation  is  running,  the  user  has  the  ability  to  interact  with  the 
model  through  the  use  of  the  following  keyboard  commands: 

•  <  P>  -  pressing  P  for  "Pause"  will  temporarily  freeze  the  simulation,  allowing  the 
user  time  to  study  the  graphical  display. 

•  <C>  -  pressing  C  for  "Continue"  counteracts  pressing  <P>  and  resumes  the 
simulation  where  it  left  off. 

•  <B>  -  pressing  B  for  "Beep"  will  turn  off  the  audible  tones  that  accompany 
customer  motion  through  the  network.  Pressing  <B>  a  second  time  will  turn 
the  beeps  back  on. 

•  <D>  -  pressing  D  for  "Display"  replaces  the  table  of  estimated  quantities  with 
the  state-versus-time  diagram  for  service  center  2.  Pressing  <D>  a  second  time 
returns  the  graphics  screen  to  the  original  configuration. 

•  <F>  -  pressing  F  for  "Fast",  accelerates  the  simulation,  causing  the  estimated 
values  to  converse  to  the  theoretical  values  more  quickly.  In  the  fast  mode,  the 
pictorial  network  display  and  state-versus-timc  graph  lor  service  center  1  are 
erased  and  not  updated.  Additionally  the  estimates  are  onlv  updated  alter  every 
10  event  occurrences.  Pressing  <  F^  a  second  time  slows  the  simulation  and 
resumes  updating  the  graphics.  Upon  return  to  the  slow  mode,  it  mav  take  some 
time  for  the  pictorial  network  display  to  catch  up  with  the  actual  status  of  the 
network. 
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•    <  S>  -  pressing  S  for  "Stop",  returns  the  simulation  to  the  program  menu. 
Because   JACKQUE    only    looks    for   keyboard   commands    while   incrementing   the 
simulation  clock,  there  may  be  a  slight  delay  between  the  user  input  and  the  program 
response. 

H.       PROGRA.M  LIMITATIONS 

JACKQUE  always  starts  the  simulation  with  the  network  completely  empty  of 
customers.  There  is  no  way  of  starting  the  simulation  in  any  other  state.  It  is  also  not 
possible  to  change  the  model  parameters  once  the  simulation  is  started.  If  it  is  desired 
to  do  so,  the  present  simulation  must  be  terminated,  the  parameters  changed,  and  a 
new  simulation  started. 

JACKQUE  uses  string  arrays  to  keep  track  of  the  customers  and  their  feedback 

count  as  they  move  through  the  network.    The  maximum  length  of  a  string  in  BASIC 

is  256.    Because  of  this  restriction,  JACKQUE  can  not  have  more  than  256  customers 

waiting  for  either  service  center  1  or  service  center  2.    If  this  condition  is  met,  the 

simulation  will  automatically  be  terminated  and  the  following  message  displayed  on  the 

screen: 

Program  execution  HALTED.  Queue 
length  exceeds  program  capabilities. 
Press  any  key  to  return  to  Main 
Program  Menu. 

JACKQUE  is  incapable  of  explicitly  tracking  any  individual  customer  for  more 
than  9  feedbacks  from  service  center  2  to  service  center  1.  Any  customer  who  has 
experienced  feedback  nine  or  more  times  will  be  represented  as  a  purple  ball  with  a 
number  "9"  in  the  middle.  Upon  seeing  this  symbol,  it  is  not  known  if  that  customer 
has  experienced  feedback  9  or  more  times.  This  limitation  causes  the  estimates  of  S^^ 
and  S2  to  be  smaller  than  the  theoretical  values.  This  bias  can  be  minimized  through 
the  judicious  selection  of  the  model  parameters  as  discussed  in  Section  IV-I. 

Although  model  events  are  occurring  continuously  in  time,  JACKQUE  updates 
the  graphics  display  only  after  the  passage  of  a  discrete  amount  of  time.  At.  It  is 
therefore  possible  for  one  or  more  events  to  occur  during  At.  When  JACKQUE  plots 
the  events  on  the  pictorial  display,  all  events  occurring  during  At  are  plotted  as  if  they 
occurred  at  the  end  of  At.  They  are  plotted  one  immediately  after  the  other,  in  the 
proper  chronological  order  in  which  they  actually  occurred. 

The  time  it  takes  the  computer  to  physically  plot  the  movement  of  the  balls  is 
not  considered  by  JACKQUE  as  part  of  the  simulation  time.    Because  the  number  of 
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events  plotted  after  each  At  is  variable,  and  because  certain  events  take  longer  to  plot 
than  others,  the  illusion  of  real  time  passage  to  the  user  is  some  what  distorted.  By 
distorted,  it  is  meant  that  the  increments  of  the  clock  as  shown  on  the  graphics  display 
are  not  always  of  equal  time  duration. 

The  illusion  of  real  time  passage  can  be  maintained,  and  the  frequency  of 
multiple  state  jumps  can  be  minimized  through  the  judicious  selection  of  the  model 
parameters  as  discussed  in  the  next  section. 

I.        GUIDANCE  ON  PARAMETER  SELECTION 

The  behavior  of  the  network  is  controlled  by  the  five  parameters  X,  ^  ,  \i  ,  P  , 
and  P^.  iMany  of  the  shortcomings  of  the  graphical  display  can  be  overcome  by  the 
careful  selection  of  these  parameters. 

The  illusion  of  real  time  passage  can  be  preserved  by  limiting  the  frequency  with 
which  multiple  events  occur  during  At.  The  inter-arrival  rate  and  service  rates  should 
be  selected  small  enough  so  that  the  probability  of  numerous  arrivals  or  service 
completions  occurring  in  At  is  small.  Picking  small  rates  will  also  improve  the  users 
ability  to  follow  the  progress  of  a  customer  in  the  network.  If  the  customers  are 
moving  too  quickly,  the  graphics  display  can  become  confusing  to  the  user.  Finally, 
the  selection  of  small  rates  will  minimize  the  number  of  times  that  the  state-versus-time 
graphs  will  exhibit  multiple  jumps. 

Although  some  feedback  makes  the  model  more  interesting,  it  is  desirable  to  keep 
the  average  amount  of  looping  for  a  customer  below  nine.  This  minimizes  the  amount 
of  bias  introduced  into  the  estimates  for  S^^  and  S^.  It  also  reduces  the  number  of 
times  a  customer  is  represented  as  a  ball  with  a  number  "9"  in  the  center.  This  is 
desirable  because  the  user  does  not  know  the  true  meaning  of  a  number  nine.  The 
amount  of  looping  is  completely  controlled  by  the  values  selected  for  P^^  and  P^. 

It  is  also  desirable  to  keep  the  average  queue  size  in  front  of  the  service  centers 
close  to  the  number  that  can  physically  be  shown  on  the  computer  screen.  Although 
this  is  not  essential,  it  makes  the  graphic  simulation  easier  to  understand.  This  is 
accomplished  by  keeping  the  service  rates  at  the  service  centers  some  what  larger  than 
the  eflective  arrival  rates. 

Selecting  the  parameters  to  accomplish  all  of  these  things  can  be  difficult  because 
they  are  so  highly  interrelated.  Table  1  contains  some  parameter  values  that  have 
yielded  interesting  results.  One  possible  method  for  selecting  parameters  is  presented 
below: 
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1.  Select  X,  the  arrival  rate  of  customers  arriving  at  service  center  1  from  outside 
the  network. 

2.  Select  P,  and  P  ,  the  probabilities  of  exiting  the  network  at  service  completion. 

3.  Compute  R^  and  R^,  the  effective  arrival  rates: 

R^  =  Vd-d-PiXi-Pp} 

R^  =  (I-P^XRi) 

4.  Select  n^  and  \i^  so  that: 

n,  >  R, 

5.  Run  the  simulation  with  these  parameters  and  see  how  it  performs. 

It  is  required  that  ^^^  and  n^  be  greater  than  zero   so  that  when  calculating  the 
theoretical  performance  quantities,  division  by  zero  is  not  attempted. 


TABLE  I 

SUGGESTED  PARAMETERS 

X 

Pi 

Pz 

^»l 

^z 

Si 

h 

Li 

Lz 

0.01 

0.10 

0.10 

0.06 

0.05 

5.3 

4.7 

7.14 

18.0 

0.05 

0.20 

0.20 

0.14 

0.12 

2.8 

2.2 

24.9 

12.5 

0.08 

0.10 

0.10 

0.60 

0.50 

5.3 

4.7 

2.35 

3.13 

0.10 

0.20 

0.30 

0.25 

0.20 

2.3 

1.8 

10.0 

10.0 

0.10 

0.10 

0.10 

0.60 

0.55 

5.3 

4.7 

7.14 

6.20 

0.10 

0.50 

0.50 

0.15 

0.07 

1.3 

.67 

8.00 

20.0 

0.15 

0.90 

0.90 

0.16 

0.05 

1.0 

.10 

17.9 

0.43 

0.20 

0.30 

0.20 

0.50 

0.35 

2.3 

1.6 

lO.O 

10.0 

0.30 

0.20 

0.60 

0.60 

0.36 

1.5 

1.2 

2.78 

50.0 

The  Jackson  network  being  modeled  can  be  transformed  into  a  simple  M/M/1 
queue  by  selecting  P^  =  1.  Likewise,  it  can  be  transformed  into  a  simple  M/M/2  queue 
by  selecting  P^  =  0  and  P^  =   1. 
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J.        PROGRAMING  MODIFICATIONS 

This  section  is  intended  to  describe  some  minor  changes  that  can  be  made  to  the 
program  JACKQUE.BAS.  A  BASIC  programing  manual  should  be  consulted  for  the 
specifics  of  how  to  change  and  then  save  the  altered  program. 

Program  line  3306  specifies  the  default  drive  on  which  the  program  looks  for  the 
random  number  generator  sub-routine  RNGEN.SRT.  To  change  the  default  drive, 
simply  replace  DRIVES  =  "A:"  with  the  appropriate  drive  label.  For  example,  if  it  is 
replaced  with  DRIVE$="B:",  the  program  will  automatically  look  on  drive  B  for  the 
sub-routine. 

Program  line  3455  contains  the  default  parameters.  These  are  the  parameters 
used  in  the  simulation  if  the  user  does  not  specify  something  else.  To  change  one  or 
more  of  the  default  parameters,  simply  replace  the  existing  value  with  the  desired  value. 

The  lengths  of  the  audible  tones  used  in  the  program  are  set  in  program  line 
3456.  The  variable  DURl  sets  the  duration  of  the  low  pitched  tone  corresponding  to  a 
customer  arrival.  The  variable  DUR2  sets  the  duration  of  the  high  pitched  tone 
corresponding  to  a  departure  from  the  network.  Either  tone  duration  can  be  changed 
by  simply  assigning  a  different  positive  integer  to  the  appropriate  variable.  The 
variable  assignment  must  also  be  changed  in  program  lines  8605  and  8610. 

The  actual  tones  of  the  sounds  are  determined  by  the  number  appearing  in  the 
statements  "SOUND  3000,  DUR2",  and  "SOUND  2500,DURr'.  To  change  the  tones, 
simply  change  the  number  appearing  after  SOUND.  These  statements  appear  in 
program  lines  12035,  18015,  18135,  23015,  and  23135. 

The  relationship  W  >  N  in  program  line  5088  determines  how  often  the 
estimated  quantities  in  the  displayed  table  are  updated.  They  are  updated  every  time 
W  exceeds  N.  W  is  the  count  of  how  many  events  have  occurred  since  the  last  update. 
The  value  of  N  is  set  in  lines  3470,  9505,  and  9515.  The  variable  W  is  initialized  in 
program  line  statements  3475  and  5088  and  incremented  in  line  5156.  As  presently 
written,  in  the  slow  mode  the  estimates  are  updated  after  every  event  and  in  the  fast 
mode  they  are  updated  after  every  10  events. 
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V.  CONCLUSIONS 

This  paper  presented  the  development  of  a  graphic  simulation  of  a  Jackson 
network.  The  method  of  graphic  simulation  allows  the  user  insight  into  the  dynamic 
nature  of  the  Jackson  network,  not  possible  through  other  instructional  methods. 

There  are  numerous  other  systems  that  could  be  better  understood  through  the 
use  of  graphic  simulation.  A  few  that  come  to  mind  are: 

Queues  with  finite  storage. 

Queues  with  bulk  arrivals 

Queues  with  bulk  service. 

Multi-channel  queues  with  one  waiting  line. 

Multi-channel  queues  with  multiple  waiting  lines,  with  or  without  jockeying. 

Priority  queuing. 

Queues  with  balking  customers. 

It  is  hoped  that  this  effort  has  been  able  to  assist  the  reader  in  better 
understanding  Jackson  networks. 
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APPENDIX  A 
JACKQUE.BAS  PROGRAM  LISTING 


50 

GOSUB 

1000 

60 

GOSUB 

2000 

70 

GOSUB 

2100 

80 

GOSUB 

5200 

90 

GOSUB 

3300 

100 

GOSUB 

3400 

110 

GOSUB 

3500 

120 

GOSUB 

3600 

130 

GOSUB 

3900 

140 

GOTO 

5000 

10  '++++++++++++++++++++++++*++++*+++++++++++++++++++++++++++++++++++++++++*++ 
20  ■+++++♦♦++♦♦+++♦   MAIN  PROGRAM  MODULE   /  MPM   ♦++♦++++++++++♦+♦+♦++++♦+♦+++ 
30  '++++++*♦++++++++++++++++++♦++++++♦+++++++++++++++++++++♦♦+++++++++++++++++ 
35  'JACKQUE.BAS,  Ver  1.0,  G.  F.  Greene,   Aug  1986 
40   CLEAR,, 7300:  KEY  OFF:  CLS 

'check  for  basica  and  color  graphics  card 

'Title  screen 

'One  time  initialization  and  Housekeeping 

'load  user  defined  functions 

'load  random  number  generator 

'initialization  for  individual  simulation  run 

'print  main  program  menu 

'compute  theoretical  statistics 

'print  main  display 

'transfer  control  to  Clock  Module 
150  SCREEN  0:  NIDTH  80:  PRINT  "Program  Ended  at  Your  Request  ..." 

160  END 

161  ' 

1000  '=======  BASICA  S  COLOR/GRAPHICS  ADAPTER  CHECK  SBRT  /  MPM  =============== 

1005  DEF  SEG  =  0:  IF  (PEEK(&H410)  AND  &H30  I  <>  aH30  THEN  DEF  SEG:  GOTO  1030 

1010  LOCATE  3,1:  PRINT  "Sorry " 

1015  PRINT  "You  do  not  have  the  color/graphics  monitor  adapter!" 

1020  PRINT  "This  simulation  uses  graphics  and  requires  that  adapter." 

1025  DEF  SEG:  END 

1030  ON  ERROR  GOTO  1035:  PLAY  "pl6":  GOTO  1050 

1035  WIDTH  80:  CLS:  LOCATE  3,1 

1040  PRINT  "This  simulation  uses  advanced  BASIC." 

1045  PRINT  "Return  to  DOS  by  typing  SYSTEM'  and  then" 

1046  PRINT  "reload  this  program  after  using  the  command  BASICA' .":END 

1050  ON  ERROR  GOTO  0:  RETURN    'return  to  line  60 

1051  ' 

2000  '====================  TITLE  SCREEN  /  MPM  ================================= 

2005  SCREEN  0,1:  COLOR  14,1,1:  WIDTH  40:  CLS 

2010  LOCATE  2,12:  PRINT  "GRAPHIC  SIMULATION" 

2015  LOCATE  4,19:  PRINT  "of  a" 

2020  LOCATE  6,12:  PRINT  "JACKSON  QUEUE  NETWORK" 

2025  LOCATE  9,14:  PRINT  "by  G.  F.  Greene" 

2030  LOCATE  14,5:  PRINT  "Submitted  in  partial  fulfillment" 

2035  LOCATE  15,2:  PRINT  "of  the  requirements  for  the  degree  of" 

2040  LOCATE  17,1:  PRINT  "Master  of  Science  in  Operations  Research" 

2045  LOCATE  19,4:  PRINT  "from  the  Naval  Postgraduate  School" 

2050  LOCATE  20,11:  PRINT  "Monterey,  California" 

2055  LOCATE  22,4:  PRINT  "Advisors  J.D.  Esary,  G.F.  Lindsay" 

2060  LOCATE  25,2:  PRINT  "Press  any  key  to  continue. .."> :  GOSUB  29000:RETURN 

2061  ' 

2100  '=======   ONE  TIME  INITIALIZATION  AND  HOUSEKEEPING  SBRT  /  MPM  ============ 

2111  SCREEN  0,1:  COLOR  14,1,1:  WIDTH  40:  CLS 

2112  LOCATE  10,5:  PRINT"Program  Initialization  ..." 

2113  LOCATE  12,13:  PRINT  "Please  Wait" 

2114  LOCATE  13,13,0   'hide  the  curser 

2120  DEFINT  A-Z   'define  all  variables  as  integers 

2130  DIM  STAT!(2),  Q(2),  NSCT!(2I,  EXIT(2),  STRG1$(255),  STRG2$( 255 )  ,SERV2( 2200 ), 

CIRC0(35),  CIRC1(35),  CIRC2(35),   CIRC3(35),  CIRC4( 35  ),circ5  (35),  CIRC6(35), 

CIRC7(35),  CIRC8(35),  CIRC9( 35  )  ,TEMPCIRC( 35 ) 
2140  DIM  ARROW1(300>,  ARR0W1C( 300 ) ,  ARROW2(100),  ARROW2C(100),  ARROW3(300), 

ARROW3C(300),  ARROW4(200),  ARR0W4C( 200  ) ,  ARROWS(600),  ARR0W5C( 600  )  ,SHIFT( 2200 ) , 

STATISTICS(2200),  SERVBYKIO),  SERVBY2(10) 
2150  FMT1$  =  "ft»»»":  FMT2$  =  "»»«».»»»"  :  FMT3$  =  "».«»":  FMT4$  =  "»».»»": 

FMT5$  =  "«»»" 
2160  TRUE  =  -1  :  FALSE  =  0:  FIRSTSCREEN  =  TRUE 
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2165 
2166 
5200 
520<+ 
3205 
3210 

3220 
3221 
3300 
3305 
3306 
3310 
3315 
3320 
3325 
3327 
3329 
3330 
3335 
3336 
33^5 
3350 
3355 
3356 
3^00 
3't20 
3^25 
3430 
3440 
3445 
3446 
3447 
3448 
3449 
3450 
3451 
3452 
3453 
3454 
3455 
3456 
3457 
3458 
3459 
3460 
3461 
3462 
3463 
3464 
3465 
3469 
3470 
3475 
3485 
3486 
3500 
3505 
3510 
3515 
3520 
3525 
3535 
3540 
3545 
3547 
3550 
3555 
3560 


RETURN 


'return  to  line  80 


<==r===============  USER  DEFINED  FUNCTIONS  SBRT  /  MPM  ==================== 

DEF  FNROUNDUPX(A!  )  =  INT(A!)  +  I 

DEF  FNEXP0NT!(A1! ,  A2! )  =  -LOG( Al!  )/A2! 

DEF  FNMINKA!  ,B!,C!  )  =  (-(A!<B!  AND  A!<C!  )*A!  )  +  (-(B!<A!  AND  B!<C!  )»B!  )  ♦ 

(-(C!<A!  AND  C»<B!  )*C!  ) 
RETURN   'return  to  line  90 


•================  LOAD  RANDOM  NUMBER  GENERATOR  /  MPM  ===================== 

RNGEN=0:  U!=0;  'initialize  random  number  generator 

DRIVE$  =  "A:"  "program  looks  on  A  drive  for  RNGEN.SRT  first 

ON  ERROR  GOTO  3329  'RNGEN.SRT  cannot  be  found  on  default  disk  drive. 

DEF  SEG  =  8H1A00:  BLOAD  DRIVE$  +  "rngen.srt",0   'load  rngen.srt 

ON  ERROR  GOTO  0   'turn  error  trapping  off 

RETURN   'return  to  line  100 

•==============FILE  LOADING  ERROR  TRAPPING  SBRT  /  MPM  ================== 

IF  ERR  <>  53  AND  ERR  <>  76  THEN  ON  ERROR  GOTO  0   'turn  error  trapping  off 

CLS:  LOCATE  3,1:  PRINT  "On  which  drive  can  the  random  number" 

PRINT  "generator  program  be  found?" 

GOSUB  29000  'capture  keyboard  input 

DRIVE$  =  RESPONSE$  +  ":" 

RESUME  3315 


FOR  I  =  1  TO 
Q(I  )  =  0 
STATd  )  = 
NSCT ! (I ) 

NEXT 

FOR  I  =  1  TO 
SERVBYKI) 
SERVBY2(I) 

NEXT 

MODE  =  2 

TIME  =  0 

CTIME!  =  0 

TAl  =  0 

TA2  =  0 

LAMBDA!  =  .1 

DURl  =  3 

UNSTABLE 

El  =  0 

E2  =  0 

ET  =  0 

NAl  =  0 

FBI  =  0 

FAST  =  0 

CUMSERVIF 

STRG1$  = 

NEXTEVENT 

N  =  1 

W  =  0 

RETURN 


INITIALIZATION  OF  SINGLE  SIMULATION  RUN  SBRT  /  MPM  ======= 

2 

'queue  counters 

'status  of  servers:  l=busy  0=idle 

'set  next  service  time  at  1  equal  to  infinity 


0 
=  IE +34 


10 


'number  of  customers  served  I  times  by  server  1/2 


'graphics  mode:  2=with  stats  table)  l=w/o  table 

'integer  clock  time 

'continuous  real  time 

'total  arrivals  at  server  1 

'total  arrivals  at  server  2 
1  :  Ul!  =  .2  :  U2!  =  .13:  PI!  =.3 
DUR2  =  2       'duration  of  audible  beeps 
0  'unstable  parameters  0=no  l=yes 

'total  exits  after  server  1 

'total  exits  after  server  2 

'ET  =  El  +  E2 

'exogenous  arrivals  at  server  1 

'feed-back  arrivals  at  server  1 

'display  speed:  O=slow  l=fast 


P2!  =.5  'default  parameters 


=  0! 


=  0 


CUMSERV2!  =  0 
STRG2$  =  "" 
XX  =  23 


'amount  of  time  server  1  S  2  busy 
'queue  strings  for  server  1  and  2 

'how  frequently  stats  table  updated 
' event  counter 


'return  to  line  110 


•=================   MAIN  PROGRAM  MENU  SBRT  /  MPM  ====== 

BORDER$  =  STRING$( 40,223) 

SCREEN  0,1:  COLOR  14,3,0:  HIDTH  40:  CLS 

LOCATE  3,14:  PRINT  "PROGRAM  MENU":  PRINT  BORDER$ 

LOCATE  6,3:  PRINT  "<I>nstructions" 

LOCATE  8,3:  PRINT  "<N>ew  random  number  Generator  Seed" 

LOCATE  10,3:  PRINT  "<D>efault  Model  Parameters" 

LOCATE  12,3:  PRINT  "<C>hange  Model  Parameters" 

LOCATE  14,3:  PRINT  "<E>nd  the  Program" 

LOCATE  15,3:  PRINT  BORDERS 

LOCATE  18,3:  PRINT  "Enter  your  selection: "> 

GOSUB  29000  'capture  keyboard  response 

AX  =  INSTRC'INDCEindce",  RESPONSES) 
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3565 

3566 
3600 
3605 
3610 
3620 
3625 
3630 
3640 
3650 
3660 
3670 
3680 
3690 
3700 
3710 
3720 
3721 
3900 
3905 
3910 
3915 
3920 
3935 
39<+0 
3945 
3950 
3955 
3960 
3965 
3970 
3980 
3981 
4000 
4001 
4002 
4003 
4004 
4006 
4010 
4012 
4014 
4016 
4018 
4020 
4022 
4024 
4025 
4028 
4030 
4032 
4034 
4036 
4038 
4040 
4050 
4051 
4052 
4054 
4056 
4058 
4060 
4062 
4064 
4065 
4066 
4067 
4068 


IF  AX  =  0  THEN  3510  ELSE  ON  AX  GOTO  4000,4200,120,4500,150,4000,4200 
,120,4500,150 


'explosive  queue 


•===========r==  THEORETICAL  STATISTICS  SBRT  /  MPM  = 

IF  (P1!=0  AND  P1!=P2!)  THEN  UNSTABLE  =1:  GOTO  4600 

Rl!  =  LAMBDA!/(1-(1-P1! )»(1-P2! )) 

R2»  =  (1  -  PI! )  *  Rl! 

IF  (R1!>=U1!  OR  R2!>=U2»)  THEN  UNSTABLE=  1:  GOTO  4600  "explosive  queue 

SI!  =  1/(P1!  ♦  ( (1-Pl!  )*(P2! ))) 

S2!  =  (1-Pl!  )*  SI! 

LI!  =  R1!/(U1!  -  Rl! ) 

L2!  =  R2!/(U2!  -  R2! ) 

L!  =  LI!  +  L2! 

Wl!  =  L1!/R1! 

W2!  =  L2!/R2! 

W!  =  L!/LAMBDA! 

Z!  =  (R1!/U1! )*(1-R2!AJ2!  ) 

RETURN   'return  to  line  130 

I 

•===============  PRINT  MAIN  DISPLAY  SBRT  /  MPM  ======================== 

SCREEN  1,0:  COLOR  0,1 

IF  FIRSTSCREEN  THEN  GOSUB  4700:  FIRSTSCREEN  =  FALSE    'one  time  printing 
GOSUB  4800  'main  screen  drawing  sbrt  of  graphics  symbols 

GET  (158,93)  -  (319,199),  STATISTICS    'put  stats  table  on  display 
PRINT  USING  FMT1$>  TIME 
PRINT  USING  FMT4$t  LAMBDA! 
PRINT  USING  FMT4§>  Ul! 
PRINT  USING  FMT4$i  U2! 
3,26:  PRINT  USING  FMT3$>  PI! 
3,37:  PRINT  USING  FMT3$>  P2! 
7,27  :  PRINT  USING  FMT5$>  Q( 2 ) 
7,10  :  PRINT  USING  FMT5$>  Q( 1 ) 


LOCATE  11,21: 
LOCATE  3,8  : 
LOCATE  9,19: 
LOCATE  9,33: 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
RETURN 


'return  to  line  140 


•===================  USER  INSTRUCTION  SBRT  /  MPM  ============== 

SCREEN  0,1:  COLOR  14,1,1:  HIDTH  40:  CLS 

LOCATE  1,6:  PRINT  "Instructions  for  the  Program  Menu" 

LOCATE  3,2:  PRINT  "Press  <s>  to  set  the  random  number  seed" j 

LOCATE  4,1:  PRINT  "If  <d>  or  <c>  is  pressed  before  setting" 

LOCATE  5,1:  PRINT  "a  seed,  the  default  seed  will  be  used" 

LOCATE  7,2:  PRINT  "Press  <c>  to  change  model  parameters." 

LOCATE  8,1:  PRINT  "The  program  will  ask  for  the  arrival" 

LOCATE  9,1:  PRINT  "rate  (lambda),  the  service  rate  for" 

LOCATE  10,1: PRINT  "server  1  (ul),  the  service  rate  for" 

LOCATE  11,1:PRINT  "server  2  (u2),  the  probability  of  exit" 

LOCATE  12,1: PRINT  "after  server  1  (PI),  and  the  probability" 

LOCATE  13,1: PRINT  "of  exit  after  server  2.   The  rates  can" 

LOCATE  14,1: PRINT  "be  any  positive  non-zero  number  while  " 

LOCATE  15,1: PRINT  "PI  S  P2  must  be  between  0  and  1.  The" 

LOCATE  16,1:PRINT  "simulation  starts  after  entering  P2." 

LOCATE  18,2: PRINT  "Press  <d>  to  use  the  default  parameters" j 

LOCATE  19,1:PRINT  "The  simulation  will  immediately  start") 

LOCATE  20,1:PRINT  "with  the  values  lambda=.10,  ul=.20," 

LOCATE  21,1: PRINT  "u2=.13,  Pl=.30,  and  P2=.50." 

LOCATE  23,2:PRINT  "Press  <e>  to  end  the  program."} 

LOCATE  25,2:PRINT  "Press  any  key  to  continue..."}:  GOSUB  29000 

CLS:  LOCATE  1,13:  PRINT  "Keyboard  Commands" 

LOCATE  3,2:  PRINT  "During  the  simulation  you  may  use  the" 

LOCATE  4,2:  PRINT  "following  keys:" 

LOCATE  6,4:  PRINT  "<p>  -  Pause  the  simulation." 

LOCATE  8,4:  PRINT  "<c>  -  Continue  thie  simulation." 

LOCATE  10, 4: PRINT  "<b>  -  turn  on  and  off  the  audio  Beeps" 

LOCATE  12, 4: PRINT  "<d>  -  change  the  Display" 

LOCATE  13, 10: PRINT  "configuration  back  and  forth  " 

LOCATE  14, 10: PRINT  "between  the  state  diagrams  and" 

LOCATE  15,10:PRINT  "the  statistics  table." 

LOCATE  17,4:  PRINT  "<f>  -  turn  on  and  off  the" 

LOCATE  18,10:  PRINT  "fast  display  mode." 

LOCATE  20,4:PRINT  "<s>  -  Stop  the  simulation  and  " 
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^070 
^♦072 
4080 
4081 
4082 
4084 
4086 
4088 
4090 
4092 
4094 
4096 
4098 
4100 
4102 
4104 
4106 
4120 
4122 
4124 
4126 
4128 
4130 
4152 
4134 
4136 
4138 
4140 
4142 
4144 
4146 
4148 
4150 
4152 
4160 
4162 
4164 
4166 
4168 
4170 
4172 
4174 
4176 
4178 
4180 
4182 
4184 
4186 
4188 
4190 
4192 
4194 
4196 
4198 
4199 
4200 
4205 
4210 
4215 
4220 
4225 
4230 
4235 
4240 
4245 
4250 
4255 
4260 
4265 


LOCATE  21>10:PRINT  "re-tum  to  the  program  menu." 
LOCATE  25,2:PRirfr  "Press  any  key  to  continue. .."> :  GOSUB  29000 
PRINT  "Graphics  Display" 


CLS:  LOCATE 
LOCATE  3,2 
LOCATE  4,1 
LOCATE  6,1 
LOCATE 


1,13: 

PRINT  "The  variables  appearing  on  the  display 


PRINT  "are  defined  as:" 

PRINT  "lambda:  Arrival  rate  from  outside  the" 

7,7:  PRINT  "queue  network." 

LOCATE  9,1:  PRINT  "ul:  Service  rate  for  server  1." 

LOCATE  11,1: PRINT  "u2 :  Service  rate  for  server  2." 

LOCATE  13,1: PRINT  "PI:  Probability  of  exit  after  server  1." 

LOCATE  15,1:PRINT  "P2:  Probability  of  exit  after  server  2." 

LOCATE  17,1: PRINT  "Ql:  Number  of  people  waiting  for" 

LOCATE  18, 5: PRINT  "server  1." 

LOCATE  20,1: PRINT  "Q2:  Number  of  people  waiting  for" 

LOCATE  21, 5: PRINT  "server  2." 

LOCATE  23,1:PRINT  "Time:  Current  time  on  simulation  clock." 

"Press  any  key  to  continue..."}:  GOSUB  29000 
PRINT  "N(t):  No.  of  people  waiting  and/or  being 

PRINT  "served  by  given  server  at  time  t." 

PRINT  "Rl:  Effective  arrival  rate  at  server  1," 

PRINT  "including  feedback." 

PRINT  "R2:  Effective  arrival  rate  at  server  2," 


LOCATE  25, 2: PRINT 
CLS:  LOCATE  1,1: 
LOCATE  2,7: 
LOCATE  4,1: 
LOCATE  5,5: 
LOCATE  7,1: 


LOCATE  8,5:  PRINT  "including  feedback." 

LOCATE  10,1: PRINT  "SI:  Expected  no.  of  times  a  person  will" 

LOCATE  11, 5: PRINT  "be  serviced  by  server  1." 

LOCATE  13,1:PRINT  "S2:  Expected  no.  of  times  a  person  will" 

LOCATE  14, 5: PRINT  "be  serviced  by  server  2." 

LOCATE  16,1:PRINT  "LI:  Long  run  expected  no.  of  people" 

LOCATE  17,5:PRINT  "either  waiting  or  being  served  by" 

LOCATE  18,S:PRINT  "server  1." 

LOCATE  20,1: PRINT  "L2:  Long  run  expected  no.  of  people 

LOCATE  21, 5: PRINT  "either  waiting  or  being  served  by" 

LOCATE  22,5:PRINT  "server  2." 

LOCATE  25,2: PRINT  "Press  any  key  to  continue. .." j :  GOSUB  29000 

CLS:  LOCATE  1,1:  PRINT  "L:   L  =  LI  +  L2" 

LOCATE  3,1:  PRINT  "HI:  Expected  time  a  person  spends" 

LOCATE  4,5:  PRINT  "waiting  or  being  served  by  server  1." 

LOCATE  6,1:  PRINT  "N2:  Expected  time  a  person  spends" 

LOCATE  7,5:  PRINT  "waiting  or  being  served  by  server  2." 

LOCATE  9,1:  PRINT  "W:   Total  expected  time  in  network." 

LOCATE  11,1: PRINT  "Z:       Long  run  proportion  of  time  server  1" 

LOCATE  12, 5: PRINT  "is  busy  and  server  2  is  idle." 

LOCATE  14,1:PRINT  "Limit:  Values  listed  under  the  column" 

LOCATE  15,8:PRINT  "heading   limit'  are  the  limiting" 

LOCATE  16,8:PRINT  "values  as  time  goes  to  infinity" 

LOCATE  17,8: PRINT  "of  the  above  defined  quantities." 

LOCATE  19,1:PRINT  "Estimate:  Values  appearing  under  the" 

LOCATE    20,ll:PRINT"hieading     estimate'    are   the" 

LOCATE  21,ll:PRINT"estimates  of  the  limiting" 

LOCATE  22,ll:PRINT"values.  They  are  calculated" 

LOCATE  23,ll:PRINT"with  data  collected  from" 

LOCATE  24,ll:PRINT"the  simulation. " > 

LOCATE  25,2:PRINT  "Press  any  key  to  return  to  the  menu ...">: GOSUB  29000 

GOTO  3500    ' return  to  the  menu 

<=========  RANDOM  NUMBER  GENERATOR  SEED  SBRT  /  MPM  ====================== 

COLOR  14,0,0:  CLS 

LOCATE  1,3:  PRINT  "Set  Random  Number  Generator  Seed":  PRINT  BORDER$ 

LOCATE  4,2:  PRINT  "Permissible  seed  values  are  integers" 

PRINT  "  in  the  range:  1  to  2147483646." 

LOCATE  7,2:  INPUT  "Enter  the  seed  value  :   "  ,A» 

IF  A»<1  OR  A»>2147483646«  THEN  LOCATE  7,2:  PRINT  STRING$( 39,32 ) :  GOTO  4225 

DEF  SEG  =  SHIAOO 

A2!  =  INT( (A»-Al!»16777210»)/65536!  ) 
into  RNGEN's  seed  storage. 


Al!  =  INT(A»/16777210»): 
'poke  the  seed's  upper  2  bytes 
POKE  8H164,A1!:  POKE  SH163,A2! 
A»=A»-A1»*16777210»-A2!»65536! : 
'poke  the  seed's  lower  2  bytes 
POKE  &H162,A1!:  POKE  iH161,A2! 


Al!  =  INT(A»/256):  A2!  =  A»-A1!»256 
into  RNGEN's  seed  storage. 


60 


<t271 
4500 
4505 
4510 
4515 
4520 
4525 
4550 

4532 
4535 
4540 

4542 
4545 
4550 

4555 
4560 
4565 

4567 
4570 
4575 

4577 
4580 
4581 
4600 
4605 
4610 
4615 
4620 
4625 
4630 
4635 
4640 
4641 
4645 
4650 
4655 
4660 
4665 
4670 
4671 
4700 
4701 
4702 
4704 
4706 
4708 
4710 
4712 
4714 
4716 
4717 
4718 
4720 
4722 
4724 
4726 
4728 
4730 
4732 
4734 
4736 
4738 
4740 


GOTO  3500 


' return  to  the  menu 


•===r==r===s==  CHANGE  MODEL  PARAMETERS  SBRT  /  MPM  ====================r= 

SCREEN  0,1:  COLOR  15,5,0:  CLS 

MSG1$  =  "Enter  a  positive,  nonzero  number  only." 

MSG2$  =  "Enter  a  positive  value  between  0  and  1." 

LOCATE  1,6:  PRINT  "CHANGE  MODEL  PARAMETERS":  PRINT  BORDER* 

LOCATE  4,1:  INPUT  "Enter  Arrival  Rate  :   ",  LAMBDA! 

IF  LAMBDA!<=0  THEN  LOCATE  23,1:PRINT  MSG1$} : LOCATE  4,1:PRINT  STRING$( 40,32 ) 

:GOTO  4525 
LOCATE  23,1:  PRINT  STRING$( 40,32  )  v 

LOCATE  6,1:  INPUT  "Enter  Service  Rate  for  Server  1  : 
IF  Ul!  <=0  THEN  LOCATE  23,1:  PRINT  MSGl$j:  LOCATE  6,1: 

GOTO  4535 
LOCATE  23,1:  PRINT  STRING$( 40,32 ) I 

LOCATE  8,1:  INPUT  "Enter  Service  Rate  for  Server  2  : 
IF  U2!  <=0  THEN  LOCATE  23,1:  PRINT  MSG1$>:  LOCATE  8,1:  PRINT  STRING$( 40,32 ) : 

GOTO  4545 
LOCATE  23,1:  PRINT  STRING$( 40,32 ) } 

LOCATE  10,1:  INPUT  "Enter  Prob  of  Exit  after  Server  1:  ",P1! 
IF  ((PI!  <  0)  OR  (PI!  >  D)  THEN  LOCATE  23,1:  PRINT  MSG2$>: 

LOCATE  10,1:  PRINT  STRING$( 40 ,32  ) :  GOTO  4560 
LOCATE  23,1:  PRINT  STRING$( 40,32  )  > 

LOCATE  12,1:  INPUT  "Enter  Prob  of  Exit  after  Server  2:  ",P2! 
IF  ((P2!  <  0)  OR  (P2!  >  D)  THEN  LOCATE  23,1:  PRINT  MSG2$>: 

LOCATE  12,1:  PRINT  STRING$( 40 ,32  ) :  GOTO  4570 
LOCATE  23,1:  PRINT  STRING$( 40,32 )) 
GOTO  120    'start  simulation 


',  Ul! 
PRINT  STRING$( 40,32)! 


U2! 


/  MPM 


LOCATE  6,5: 
LOCATE  7,5: 
LOCATE  8,5: 
LOCATE  9,5: 


•r==r=========r  EXPLOSIVE  QUEUE  WARNING  SBRT 

SCREEN  0,1:  COLOR  14,3,0:  WIDTH  40:  CLS 

PRINT  "Your  parameter  selection  results" 
PRINT  "in  an  Unstable  (EXPLOSIVE)  queue." 
PRINT  "The  graphical  presentation  of  the" 
PRINT  "network  will  still  be  accurate^" 
LOCATE  10, 5: PRINT  "However  the  network  statistics" 
LOCATE  11, 5: PRINT  "will  be  INVALID,  and  will  NOT  be" 
LOCATE  12,5:PRINT  "displayed.  See  the  user's  manual" 
LOCATE  13, 5: PRINT  "for  an  explanation." 
LOCATE  15,10:PRINT  "<C>ontinue  with  Simulation" 
LOCATE  16, 10: PRINT  "<R>eturn  to  the  Program  MENU" 
LOCATE  18,5:  PRINT  "ENTER  your  Selection" 
GOSUB  29000  'Capture  Keyboard  input 
AX  =  INSTRC'CRcr",  RESPONSE*) 
IF  A/C  =  0  THEN  4605  ELSE  ON  AX  GOTO  3720,  100,  3720,  100 


>==============  DRAW  &   SAVE  CIRCLES/ARROWS  SBRT  /  MPM  ===== 

'draw  circles 

CIRCLE  (11,35),7,1:PAINT(11,3S),1,1 

Y=4  :  X=27 

FOR  1=1  TO  9 

CIRCLE  (X,35),7,2:  PAINT( X,35) ,2,2 

LOCATE  5,Y  :  PRINT  RIGHT$(STR$(  I ),1 ) 

X=X+16 

Y=Y+2 
NEXT 

'draw  arrows  colored  and  uncolored 

DRAW  "c3>  bm  0,47j  m4,47>h4>f4ig4ie4>  ml35,47j  h4j  f4»  g4" 
DRAW  "c2>  bm  0,63>  m4,63 >h4>f4jg4>e4>  ml35,63>  h4>  f4>  g4" 
DRAW  "c3v  bm  159, 47 >  m  271, 47 j  h4i  f4>  g4" 
DRAW  "c2>  bm  159, 63j  m  271, 63>  h4>  f4>  g4" 
DRAW  "c3j  bm  303, 79 >  m  15,79>  f4>  h4>  e4" 

f4  h4j  e4>  bm  111,79)  f4j  h4j  e4" 
m  15,95j  f4>  h4>  e4" 
DRAW  "cZi   bm207,95>  f4i  h4>  e4v  bm  lll,95j  f4v  h4)  e4" 
DRAW  "c3v  bm  295, 47>  m319,47v  h4v  f4>  g4" 
DRAW  "c2>  bm  295, 63>  m  319,63)  h4j  f4j  g4" 
'c3)  bm  223, 39j  m  223,11}  f4>  h4)  g4" 


DRAW  "c3>  bm  207, 79j 
DRAW  "c2j  bm  303, 95 j 


DRAW 


DRAW  "c2j  bm  239, 39j  m  239, 11>  f4)  h4>  g4" 
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<f7^1 
^742 

4746 
4748 
4750 
4752 
4754 
4756 
4758 
4760 
4761 
4762 
4764 
4766 
4768 
4770 
4772 
4774 
4776 
4778 
4780 
4782 
4784 
4785 
4800 
4801 
4802 
4804 
4806 
4808 
4810 
4812 
4814 
4816 
4818 
4820 
4822 
4824 
4826 
4828 
4830 
4832 
4834 
4836 
4838 
4840 
4841 
4842 
4844 
4846 
4848 
4850 
4852 
4854 
4856 
4858 
4860 
4862 
4864 
4866 
4868 
4870 
4872 
4874 
4876 
4878 
4880 
4882 


'save  circles  for  future  plotting 
GET  (4,28)  -  (18,42),  CIRCO 
(20,28)  -  (34,42),  CIRCl 
(36,28)  -  (50,42),  CIRC2 
(52,28)  -  (66,42),  CIRC3 
(68,28)  -  (82,42),  CIRC4 
(84,28)  -  (98,42),  CIRC5 
(100,28)  -  (114,42),  CIRC6 
(116,28)  -  (130,42),  CIRC7 


GET 
GET 
GET 
GET 
GET 
GET 
GET 
GET 
GET 


( 132,28) 
(148,28) 


(146,42), 
( 162,42), 


CIRCS 
CIRC9 


'save  arrows  for  future  plotting 


GET(0,47)  - 
GET(0,63)  - 
GET(223,11) 
GET(239,11) 
GET( 159,47) 
GET( 159,63) 
GET( 295,47) 
GET( 295,63) 
GET(15,79)  - 
GET( 15,95)  - 
CLS 
RETURN   ' return 


(135,46), 
(135,62), 


ARROWl 
ARROWIC 


(224,39), 
(240,39), 
(271,48), 
(271,64), 
(319,48), 
(319,64), 
(303,80), 


ARR0W2 
ARR0W2C 
ARR0W3 
ARR0W3C 
ARR0W4 
ARR0H4C 
ARR0W5 


(303,96),  ARR0W5C 


to  line  3915 


<=================  MAIN  SCREEN  DRAWING  SBRT  /  MPM  ===================== 

'print  network  picture 

LINE  (0,8)-(319,8),2  :  LINE  (0,  10 )-( 319,10  ),2 

LOCATE  1,14:  PRINT  "JACKSON  QUEUE", 

DRAW  "c3j  bm  0,43>  m  135,43  J  h4j  f4>  g4" 

LINE  (135,31)-(159,55),3,B 

DRAW  "c3>  bm  159, 43j  m  271, 43>  h4»  f4>  g4" 

DRAW  "c3j  bm  167,43 >  m  167, 15 >  f4>  h4j  g4" 

LINE  (271,31)  -  (295, 55), 3, B 

DRAW  "bm  295, 43>  m  319, 43v  h4j  f4j  g4" 

DRAW  "bm  303, 43 >  m  303,76}  h4j  f4}  e4>  g4" 

DRAW  "m  15,76>  f4>  h4>  e4}  g4" 

DRAW  "m  15,43}  f4}  h4}  g4" 

DRAW  "bm4,43}  h4j  f4j  g4}  bm  207, 76>  f4}  h4)  e4>  bm  111,76}  f4}  h4}  e4" 

LOCATE  3,1:  PRINT  "lambda=" 

LOCATE  3,23:  PRINT  "Pl=" 

LOCATE  3,34:  PRINT  "P2=" 

LOCATE  7,7:  PRINT  "Ql=" 

LOCATE  9,16:  PRINT  "ul=" 

LOCATE  7,24:  PRINT  "Q2=" 

LOCATE  9,30:  PRINT  "u2=" 

LOCATE  11,16:  PRINT  "TIME:" 

'print  state-verses  time  graphs 

LINE  (0,90)  -  (319, 90), 2  :  LINE  (0,92)  -  (519, 92), 2 

LINE  (157,92)-(157,199),2 

LINE  ( 23,185)-(23,104) 

LINE  ( 183, 185  )-( 183,104) 

FOR  Y  =  107  TO  179  STEP  8 

LINE  (19,  Y)  -  (23,Y) 

LINE  (179,Y)-(183,Y) 
NEXT 

LINE  ( 23, 185)-( 151,185) 
LINE  ( 183,185  )-( 311,185) 
FOR  X  =  27  TO  151  STEP  4 

IF  (X=63)  OR  (X=103)  OR  (X=145)  THEN  GOTO  4870 

LINE  (X,185)  -  (X,188) 

LINE  (X+160,  185)-(X+160,188):  GOTO  4874 
LINE  (X,185)-(X,190) 
LINE  (X+160, 18S)-(X+160, 190) 
NEXT 
FOR  COL  =  1  TO  21  STEP  20:  K  =  10 

FOR  ROW  =  15  TO  23  STEP  2 
LOCATE  ROW,  COL 
K=K-2 


62 


4884 
4886 
4888 
4890 
4892 
4894 
4896 
4898 
4900 
4902 
4904 
4906 
4908 
4910 
4912 
4914 
4916 
4918 
4920 
4922 
4924 
4926 
4928 
4930 
4932 
4934 
4936 
4937 
5000 
5010 
5030 
5031 
5032 
5033 
5034 
5050 
5080 
5081 
5085 
5088 
5090 
5095 
5100 

5110 
5115 
5120 
5125 
5126 
5127 
5130 
5135 
5140 
5144 
5145 
5146 
5147 
5148 
5149 
5150 
5151 
5152 
5155 
5156 
5160 
5161 
5162 
6000 
6010 


PRINT  STR$(K) 
NEXT 
NEXT 

LOCATE  13,1:  PRINT  "N( t  )" 
LOCATE  13,21:  PRINT  "N( t  )" 
LOCATE  25,8:  PRINT  "SERVER  l"j 
LOCATE  25,27:  PRINT  "SERVER  2"> 
GET  (158,93)-(319,199),  SERV2 
LINE  (158,93)-(319,199),0,BF 
IF  UNSTABLE  =  0  THEN  GOTO  4910 
LOCATE  18,  24:  PRINT  "Unstable  Queue" 
LOCATE  19,23:  PRINT  "Model  Statistics" 
LOCATE  20,25:  PRINT  "Are  Invalid":  GOTO  4936 

PRINT  "VAR  ESTIMATE 

PRINT  "Rl    0.000" 

PRINT  "R2 
SI 


'save  server  2  state  graph 
'erase  server  2  state  graph 
'stable  do  not  print  warning 


LOCATE  13,21: 
LOCATE  14,21: 


LIMIT' 


0.000" 
PRINT  "SI  0.000" 
PRINT  "S2  0.000" 
PRINT  "LI  0.000" 
PRINT  "L2  0.000" 
PRINT  "L  0.000" 
PRINT  "Wl  0.000" 
PRINT  "W2  0.000" 
PRINT  "W  0.000" 
PRINT  "Z     0.000"^ 

Y  =  32:  GOSUB  30000 

RETURN    'return  to  line  3920 


LOCATE  15,21: 

LOCATE  16,21: 

LOCATE  17,21: 

LOCATE  18,21: 

LOCATE  19,21: 

LOCATE  20,21: 

LOCATE  21,21: 

LOCATE  22,21: 

LOCATE  23,21: 

LOCATE  24,21: 


'print  stat  limits  to  screen 


'+++++++♦++++++♦++++++♦++♦♦+++++++♦+♦++♦++++++++♦+♦♦♦♦+♦♦+♦++++♦+♦+♦+♦ 
'++++♦++♦++++++++♦♦+    CLOCK  MODULE   /  CM   +++++♦+++♦♦++++++++♦♦+♦+++♦+♦ 
'++++♦++++++♦+++++++++++++++++++++++♦+++++++++++♦♦+♦♦++++++++++++++++♦ 
DEF  SEG  =  SHIAOO:  CALL  RNGEN(U!)    'generate  a  uniform(0,l)  RV 
NAT!  =  FNEXPONT»(U!,  LAMBDA!)        'initialize  first  exogenous  arrival 
GOSUB  25000     'initialize  exit(l)   -  exit  after  server  1?:  l=yes  0=no 
GOSUB  26000     'initialize  exit(2)   -  exit  after  server  2?:  l=yes  0=no 
LOCATE  7,27:  PRINT  USING  FMT5$>  Q( 2 ) 
LOCATE  7,10:  PRINT  USING  FMT5$J  Q(l ) 

IF  l(Q(l)  >=  250)  OR  (Q(2)  >=  250))  THEN  GOSUB  9600   'program  crash 
IF  UNSTABLE  =  1  THEN  GOTO  5090   'skip  updating  statistics 

IF  (H>=N  AND  MODE  =  2)  THEN  H=0:  GOSUB  9000   'update  and  print  statistics 
CTIME!  =  FNMIN!(NAT!,  NSCT!(1),  NSCT!(2))  'select  next  event  time 
IF  CTIME!  =  0  THEN  GOSUB  7500    'tie-breaking  sbrt 
IF  CTIME!  =  NAT!  THEN  NEXTEVNT  =  1  ELSE  IF  CTIME!  =  NSCT!(1)  THEN 
NEXTEVNT  =  2  ELSE  NEXTEVNT  =  3    'determine  what  next  event  is 
NEXTEVNTTIME  =  FNROUNDUPXC CTIME !  )    'round  up  to  next  integer 
DEF  SEG:  POKE  8H6A,0:  'clear  keyboard  buffer 

WHILE  TIME  <  NEXTEVNTTIME 
TIME  =  TIME  ♦  1 

IF  FAST  =  1  THEN  GOTO  5145     'skip  updating  state  graphs 
'if  time  <=  32  skip  screen  shift  sbrt>  goto  state  diagram  update  sbrt 
IF  TIME  <=  32  THEN  GOSUB  6000:  XX=XX+4:  FOR  I  =  1  TO  600:  NEXT:  GOTO  5145 
GOSUB  7000   'screen  shift  subroutine  if  time  >=  32 
XX  =  147:  GOSUB  6000   'state  diagram  update  sbrt 
LOCATE  11,21:  PRINT  USING  FMT1$>  TIME 
'check  keyboard  buffer  for  user  input  and  act  accordingly 

'if  no  keyboard  input  skip 


A$  =  INKEY$:  IF  A$  =  ""  THEN  GOTO  5152 

IF  (A$  =  "s"  OR  A$  =  "S")  THEN  GOTO  100 

"p"  OR  A$  =  "P")  THEN  GOSUB  8000 

"d"  OR  A$  =  "D")  THEN  GOSUB  8500 

"B"  OR  A$="b"  )  THEN  GOSUB  8600  : 

"F"  OR  A$  =  "f")  THEN  GOSUB  9500 
POKE  &H6A,0:    'clear  keyboard  buffer 


IF  (A$  =  "p"  OR  A$  =  "P")  THEN  GOSUB  8000  :  GOTO  5152 

IF  (A$  =  "d"  OR  A$  =  "D")  THEN  GOSUB  8500  :  GOTO  5152 

IF  (A$  =  "B"  OR  A$="b")  THEN  GOSUB  8600  :  GOTO  5152 

IF  (A$  = 

DEF  SEG: 

WEND 

W  =  W+1  'event  counter  -  determines  when  stats  updated  in  line  5088 

ON  NEXTEVNT  GOTO  10000,  15000,  20000    'arrival  module,  server  1  module  or 

server  2  module 
I 

•  =========  SERVER  1  STATE  DIAGRAM  UPDATE   SBRT  /  CM   =  ================= 

CI  =  2   'determines  color  of  tic  marks:  l=blue  2=purple 
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6020 

6030 
60^0 
60^5 
6050 
6051 
6500 
6510 
6520 
6530 
65^+0 
6550 
6551 
7000 
7010 
7020 
7025 
7030 
70^0 
7050 
7051 
7500 
75'*0 
7550 
7560 
7570 
7580 
7585 
7590 
7595 
7597 
7600 
7605 
7610 
7615 
7620 
8000 
8005 
8010 
8020 
8021 
8500 
8505 
8510 
8520 
8530 
8531 
8600 
8601 
8605 
8610 
8615 
8616 
9000 
9005 
9010 
9015 
9020 
9021 
9022 
9023 
9024 
9025 
9030 
9035 
9040 
9045 
9050 
9055 


NTl  =  STAT(l)  -f  Q(  1 )    'ioial  customers  waiting  and  being  served  3  1 
IF  NTl  >  9  THEN  NTl  =  9:  CI  =  1   'can  not  plot  more  than  9  circles 
LINE  (XX,  179  -  (8*NT1))  -  STEP(4,0),  CI   'plot  the  tic  mark 
IF  M0DE=1  THEN  GOSUB  6500  'if  server  2  state  graph  is  displayed  update 
RETURN    "return  to  line  5130  or  5145 


•==r======  SERVER  2  STATE  DIAGRAM  UPDATE   SBRT  /  CM   =============== 

C2  =  2 

NT2  =  STAT(2)  *   QIZ)         'total  customers  waiting  and  being  served  3 
IF  NT2  >  9  THEN  NT2  =  9:  C2  =  1   'can  not  plot  more  than  9  circles 
LINE  (XX  ♦  160,  179  -  (8»NT2))  -  STEP(4,0),  C2   'plot  the  tic  mark 
RETURN    'return  to  line  5130  or  5145 


SCREEN  SHIFT  SBRT  /  CM 


GET  (28,  103)  -  (155,183),  SHIFT 

PUT  (24,103),  SHIFT,  PSET 

IF  MODE  =  2  THEN  RETURN 

GET  (188,103)  -  (315,183),  SHIFT 

PUT  (184,103),  SHIFT,  PSET 

RETURN     'return  to  line  5140 


'capture  server  1  state  graph 
'shift  and  re-plot  the  graph 
'if  not  displaying  server  2  return 
'capture  server  2  state  graph 
'shift  and  replot  the  graph 


DEF  SEG 
IF  NSCT! 
IF  NAT! 
IF  NAT! 
IF  U!  <= 
CTIME!  = 
IF  U!  <= 
IF  U!  >= 
CTIME!  = 
IF  U!  <= 
CTIME!  = 
IF  U!  <= 
CTIME!  = 


==========   TIE  BREAKING  SBRT  /  CM  ==== 

=  SHIAOO:  CALL  RNGEN(U!) 
(1)  =  NSCT!(2)  AND  NSCT'd)  =  NAT!  GOTO 
=  NSCT!(2)  GOTO  7600       '2-viay  tie 
=  NSCT!(1)  GOTO  7610       '2-way  tie 
.5  THEN  CTIME!  =  NSCT!(1):  NEXTEVNT 
NSCT!(2):  NEXTEVNT  =  3:  GOTO  5110 
.333  THEN  CTIME!  =  NAT!:  NEXTEVNT  = 
.333  AND  U!  <=  .666  THEN  CTIME!  = 
NSCT!(2):  NEXTEVNT  =  3:  GOTO  5110 
.5  THEN  CTIME!  =  NAT!:  NEXTEVNT  = 
NSCT!(2):  NEXTEVNT  =  3 
.5  THEN  CTIME!  =  NAT!: 
NSCT!(1):  NEXTEVNT  =  2 


7590 


'3-way  tie 


:  GOTO  5110 

NEXTEVNT  = 

:  GOTO  5110 


=  2:  GOTO  5110  'return  clock 

=  1:  GOTO  5110  'return  to  clock 
NSCT!(1):NEXTEVNT=2:  GOTO  5110 

' return  to  clock  module 
1:  GOTO  5110   'return  to  clock 

' return  to  clock  module 
1:  GOTO  5110   'return  to  clock 

' return  to  clock  module 


•=====================  PAUSE  SBRT  /  CM  ===================== 

'used  in  flashing  of  arrows 

DEF  SEG:  POKE  8H6A,0: 

A$  =  INKEY$:  IF  A$="c"  OR  A$="C"  THEN  RETURN  ELSE  GOTO  8020 


•=========  CHANGE  DISPLAY  CONFIGURATION  SBRT  /  CM  =============== 

'replaces  server  2  state  diagram  with  stats  table  or  vice/versa 
IF  MODE  =  1  THEN  PUT( 158,93),  STATISTICS,  PSET:  MODE  =  2:  RETURN 
PUT  (158,93),  SERV2,  PSET:  MODE  =  1 
RETURN   'return  to  line  5149 


•================   BEEP  SUPPRESSION  SBRT  /  CM 

'turns  the  audible  beeps  on  and  off 
IF  DURl  =  3  THEN  DURl  =  0  ELSE  DURl  =  3 
IF  DUR2  =  2  THEN  DUR2  =  0  ELSE  DUR2  =  2 
RETURN   "return  to  line  5150 


<=========  UPDATE  STATISTICAL  ESTIMATES  SBRT  /  CM  == 

TAl  =  NAl  +  FBI 

Rl!  =  TAl  /  CTIME! 

R2!  =  TA2  /  CTIME! 

IF  CUMSERVl!  =  0  THEN  SRI!  =  0:  SR2!  =  0:  GOTO  9024 


SRI!  =  (El  +  TA2)/CUMSERV1! 
IF  CUMSERV2!  =  0  THEN  SR2  =  0: 
SR2!  =  (E2  ♦  FB1)/CUMSERV2! 
R!  =  NAl/CTIME! 


'estimate  of  service  rate  for  server  1 
GOTO  9024 

"estimate  of  service  rate  for  server  2 
'estimate  of  exogenous  arrival  rate 


S1!=0:  S2!=0:  ET 
FOR  I  =  1  TO  10 
SI!  =  SI!  + 
S2!  =  S2!  ♦ 
NEXT 

SI!  =  S1!/ET 
S2!  =  S2!/ET 


El  ♦  E2:  IF  ET  =  0  THEN  GOTO  9056 

I»SERVBY1(I) 
I»SERVBY2(I ) 
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9056  IF  (SRI!  <=  Rl! )  THEN  LI!  =  0:  GOTO  9062 

9060  H!  =  R1!/(SR1!  -  Rl!  ) 

9062  IF  (SR2!  <=  R2! )  THEN  L2!  =  0:  GOTO  9070 

9065  L2!  =  R2!/(SR2!  -  R2! ) 

9070  L!  =  LI!  +  L2! 

9075  Wl!  =  L1!/R1! 

9077  IF  (R2!  =0)  THEN  W2!=0:  GOTO  9085 

9080  M2!  =  L2!/R2! 

9085  W!  =  L!/R! 

9086  IF  SRI!  =  0  THEN  Z!  =  0:  GOTO  9200 

9087  IF  (PI!  =  1)  THEN  Z!  =  (R1!/SR1!):  GOTO  9200 

9088  IF  SR2!  =  0  THEN  Z!  =  0:  GOTO  9200 
9090  Z!  =  ABS( (R1!/SR1! )*(1  -  (R2!/SR2!))) 

9200  Y=23:  GOSUB  30000     'plot  estimates  to  screen 

9210  RETURN    'return  to  line  5090 

9211  ' 

9500  '==============  GRAPHICS  SUPPRESSION  SBRT  /  CM  ========================== 

9501  'if  fast  mode  is  selected  no  graphics  are  displayed 
9505  IF  FAST  =  0  THEN  FAST  =  1  ELSE  FAST  =0:  N  =  1 
9510  IF  FAST  =  1  THEN  GOTO  9515  ELSE  RETURN 

9515  N  =  10  'when  fast  mode  selected)  stats  updated  after  every  N  events 

9520  LINE  (l'^0,26)  -  ( 15<+,50 )  ,0 ,BF    'erase  customers  and  state  tic  marks 

9525  LINE  (4,28)  -  ( 150 ,42  )  ,0 ,BF 

9530  LINE  (276,36)  -  ( 290 ,50  )  ,0 ,BF 

9535  LINE  (172,28)  -  ( 266 ,42  )  ,0 ,BF 

9540  LINE  (24,105)  -  ( 151 ,183  )  ,0 ,BF 

9545  RETURN    'return  to  line  5152 

9600  '================  PROGRAM  CRASH  SBRT  /  CM  ============================== 

9605  SCREEN  0,1:  COLOR  14,3,0:  WIDTH  40:  CLS 

9610  LOCATE  7,3:  PRINT  "Program  execution  HALTED.  Queue" 

9615  LOCATE  8,3:  PRINT  "length  exceeds  program  capabilities." 

9620  LOCATE  9,3:  PRINT  "Press  any  key  to  return  to  Main" 

9621  LOCATE  10,3:  PRINT  "Program  menu."j 

9625  GOSUB  29000:  GOTO  100   'display  program  menu 

9626  ' 

10000  •++++++++++++t+++++++*+++++++++++++++++++++++++++++++++++++++++++++++++ 

10010  '++++++++++++++♦+++    ARRIVAL  MODULE  /  AM    +♦++♦♦++++++++++++++♦++++++ 

10020    '+  +  +  +  +  +  +  +  +  +  +  +  ♦  +  +•  +  +  +  +  ■♦•++  +  +  +  +++  +  +  +  +  +•  +  +++  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  ♦  +  + 

10030  STRG1$  =  "0"  +  STRG1$    '0  corresponds  to  new  arrival,  add  to  queue  strg 

10035  NAl  =  NAl  +  1 

10037  IF  FAST  =  1  THEN  GOTO  10042  'skip  plotting  routine 

10040  IF  Q(l)<8  THEN  GOSUB  11000   ELSE  GOSUB  12000:  PUT( 0,42 ) ,ARR0W1 ,PSET 

10042  DEF  SEG  =  SHIAOO:  CALL  RNGEN(U!) 

10045  AINC!  =  FNEXPONT!(U! ,  LAMBDA!)  'compute  next  inter-arrival  time 

10050  NAT!  =  AINC!  +  CTIME!  'compute  new  exogenous  arrival  time 

10060  IF  STAT(l)  =  1  THEN  Q( 1 )  =  Q( 1)  +  1:  GOTO  5080    'return  to  clock  module 

10063  CALL  RNGEN(U!  ) 

10065  SlINC!  =  FNEXPONT!(U! ,U1!  )     'compute  new  server  1  service  duration 

10070   NSCT!(1)  =  SlINC!   ♦  CTIME!   'compute  next  service  completion  time  for  1 

10080  STAT(l)  =  1   'change  status  of  server  1  to  busy 

10090  GOTO  5081     'return  to  the  clock  module 

10091  ■ 

11000  ■===============   ARRIVAL  PLOTTING  SBRT  /  AM  ========================== 

11005  GOSUB  12000   'flash  arrow  1 

11010  IF  STAT(l)  =  0  THEN  PUT( 140,36),  CIRC0:GOSUB  28000:  PUT(0,42),  ARROWl,PSETi 

RETURN   'return  to  line  10040 
11020  PUT(116-Q(1)*16,28),  CIRCO: GOSUB  28000:  PUT  (0,42),  ARROWl,  PSET 

11030  RETURN     'return  to  line  10040 

11031  • 

12000  '================   FLASH  ARROW  1  SBRT  /  AM  ============================ 

12010  PUT (0,42),  ARROWIC,  PSET     'flash  purple  arrow 

12020  FOR  I  =  1  TO  500 

12030  NEXT 

12035  SOUND  500,DUR1 

12040  RETURN  'return  to  line  10040  or  11010 

12041  ' 

15000  '++♦++++++♦+++++++++++++++++++++♦+++++++++++++++++♦+♦♦+++++++♦++•+♦+++♦♦ 
15010  '++++++++++++++++  SERVICE  CENTER  1  MODULE  /  SCIM  ++++++++♦++♦+++++++++ 
15020  '+++++++♦+++++++++++♦♦++♦+♦♦++♦++♦+++++♦+♦++++++♦♦++++++++++++++++♦++++ 
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1502S 
15030 
15035 
150^0 
15050 
15060 
15070 
15072 
15075 
15080 
15090 
15099 
15100 
15110 
15113 
15115 
15120 
15130 
15140 
151<il 
17000 
17010 
17020 
17030 
17035 
170'40 
17050 
17055 
17060 
17070 

17080 
17090 
17095 
17099 
17100 
17110 
17115 
17120 
17130 
171'^0 
17150 
1715'+ 
17155 
17160 
17170 
17175 
17180 
17181 
18000 
18010 
18015 
18020 
18030 
180'+0 
18050 
18051 
18100 
18110 
18120 
18130 
18135 
181'+0 
18141 
18500 
18501 
18510 
18520 
18530 


CUMSERVl!  =  CUMSERVl!  ♦  SlINC!  'update  amount  time  server  2  busy 

TRANSFER1$  =  RIGHT$( STRG1$ ,1 )   'move  customer  from  Ql  to  Q2 

IF  EXIT(l)  =  0  THEN  STRG2S  =  TRANSFER1$  +  STRG2$:  TA2  =  TA2+1 

STRG1$  =  LEFT$(  STRG1$,  LEN(STRG1$)  -  1) 

IF  FAST  =  0  THEN  GOSUB  17000   'server  1  plotting  sbrt 

IF  EXIT(l)  =  1  THEN  El  =  El  +  1:  GOSUB  18500:  GOTO  15100   'estimate  31 

IF  STAT(2)  =  1  THEN  Q(2)  =  Q(2)  ♦  1  :  GOTO  15100 

DEF  SEG  =  SHIAOO:  CALL  RNGEN(U!) 

S2INC!  =  FNEXPONT!(U! ,  U2!  )   'compute  new  server  2  service  duration 

NSCT!(2)  =  S2INC!  +  CTIME!    'compute  new  service  completion  time 

STAT(2)  =  1    'change  status  for  server  2  to  busy 
'make   server  1  completion  impossible  event 
IF  Q(l)  =  0  THEN  NSCT!(1)  =  lE+31:  STAT(l)  =  0:  GOTO  15130 

Q(  1)  =  Q(  1)  -  1 
DEF  SEG  =  SHIAOO:  CALL  RNGEN(U!) 

SlINC!  =  FNEXPONT!(U! >  Ul!  )    "compute  new  server  1  service  duration 

NSCT!(1)  =  SlINC!  +  CTIME!     'compute  new  service  completion  time 
GOSUB  25000        'generate  new  exit( 1 )  indicator  variable 
GOTO  5050         'return  to  clock  module 

•=====rr=========   SERVER  1  PLOTTING  SBRT  /  SCIM  ====================== 

GET (140,36)  -  (154,50),  TEMPCIRC      'capture  circle  3  server  1 

LINE(140,36)  -  (154,50),0,BF         'erase  circle  3  server  1 

IF  EXIT(l)  =  1  THEN  GOSUB  18000:  GOTO  17060    'flash  arrow  2 

GOSUB  18100   'flash  arrow  3 

IF  STAT(2)  =  0  THEN  PUT( 276,36),  TEMPCIRC  :  GOTO  17055  'put  at  server  2 

IF  Q(2)<  6  THEN  PUT( 252-( Q( 2  )*16  )  ,28  ) ,  TEMPCIRC    'put  in  server  2  queue 

GOSUB  28000  : PUT ( 159,43  ) ,  ARR0W3,  PSET    'plot  white  arrow  3 

IF  Q(l)  =  0  THEN  RETURN    'return  to  line  15060 

GET(116,28)-(130,42),  TEMPCIRC  :  LINE  (116 ,28  )-(130,42 ) ,0,BF: 
PUT( 140,36),  TEMPCIRC   'put  next  customer  in  server  1  box 
IF  Q(l)  =  1  THEN  RETURN     'return  to  line  15060 
X  =  100 

IF  Q(l)<=  8  THEN  K  =  (Q(l)  -1)  ELSE  K  =  7 
'shift  all  remaining  customers  to  the  right  one 
FOR  J  =  1  TO  K 

GET(X,28)  -  (X+14,42),  TEMPCIRC 

LINE  (X,28)  -(X+14,42),0,BF 

PUT(X+16,28),  TEMPCIRC 

X  =  X-16 
NEXT 

IF  Q(l)<9  THEN  RETURN   'return  to  line  15060 

'go  to  overflow  queue  and  bring  next  customer  into  visible  queue 
L00PC0UNT1$  =  MID$(STRG1$,(LEN(STRG1$)-8),1) 
X=4  :  Y=28 

INDEX  =  VAL(L00PC0UNT1$)  +1 

GOSUB  27000   'bring  circle  in  from  queue  overflow 
RETURN   'return  to  Ine  15060 


•=================   FLASH  ARROW  2  SBRT  /  SCIM  ====== 

PUT  (167,15),  ARR0W2C,  PSET    'flash  purple  arrow  2 

SOUND  3000, DUR2 

FOR  1=1  TO  1000 

NEXT 

PUT (167,15),  ARR0W2,  PSET      'plot  white  arrow  2 

RETURN     'return  to  line  17030 


•=================   FLASH  ARROW  3  SBRT  /  SCIM  ========================= 

PUT  (159,43),  ARR0W3C,  PSET    'flash  purple  arrow  3 

FOR  1=1  TO  500 

NEXT 

SOUND  500,DUR1 

RETURN       'return  to  line  17030 

•=========  FREQUENCY  OF  SERVICE  BY  SERVER  1   SBRT  /  SCIM  =============== 

'this  sbrt  used  in  calculation  of  SI 

FOR  I  =  0  TO  9 

IF  VAL(TRANSFER1$)  =  I  THEN  SERVBYKI  +  l)  =  SERVBYKI+l)  +  1:  GOTO  18540 

NEXT 
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18535  IF  VAL(TRANSFER1$)  =  0  THEN  RETURN 

185«»0  FOR  I  =  1  TO  9 

18550  IF  VAL(TRANSFER1$)  =  I  THEN  SERVBYZd)  =  SERVBY2(I)  ♦  1:  RETURN 

18560  NEXT 

18565  RETURN      'return  to  line  15060 

18566  ' 

20000  '++•  +++++♦++♦♦  +  ♦♦  +  +  +■++++■++  +  +  +  +  ♦++++♦++++♦♦♦+++♦++  +  +  ♦  ++++♦  +  +  +  +  ♦+♦++♦+  +  +  + 

20010  •++++♦+++♦+++++  SERVICE  CENTER  2  MODULE  /  SC2M  ++♦+++++++++♦++♦++++++++ 

20020  '+  +  +  +  *  +  +  +  -t-  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  *  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  -»-  +  +  + 

20025  CUMSERV2!  =  CUMSERV2!  +  S2INC!  'update  amount  of  time  server  2  busy 

20027  TRANSFER2$  =  RIGHT$( STRG2$ ,1 ) 

20030  IF  EXIT(2)  =  1  THEN  E2  =  E2+1:  GOSUB  23500:  GOTO  20070   'estimate  S2 

200<i5   FBI    =    FBI    ■♦•    1      'increment    feedback  at   server   1   counter 

20050  IF  (VAL(TRANSFER2$)<9)  THEN  TRANSFER2$  =RIGHT$(  STR$(  VAL( TRANSFER2$ )+l ) ,1 ) 

ELSE   TRANSFER2$  =  "9" 
20060  STRG1$  =  TRANSFER2$  +  STRG1$ 
20070  STRG2$  =  LEFT$( STRG2$ ,  LEN(STRG2$)  -  1) 
20080  IF  FAST  =  0  THEN  GOSUB  22000  'server  2  plotting  subrt 
20090  IF  EXIT(2)  =  1  THEN  GOTO  20130 

20100  IF  STAT(l)  =  1  THEN  Q( 1 )  =  Q(l )  +1  :  GOTO  20130 
20102  DEF  SEG  =  8H1A00:  CALL  RNGEN(U!) 

20105  SlINC!  =  FNEXPONT!(U! ,  Ul!  )     'compute  new  server  1  service  duration 
20110  NSCT!(1)  =  SlINCJ  +  CTIME!      'compute  new  service  completion  time 
20120  STAT(l)  =  1   'change  status  of  server  1  to  busy 

20129  'make  service  completion  at  server  2  impossible  event 

20130  IF  Q(2)  =  0  THEN  NSCT!(2)  =  lE+31:  STAT(2)  =  0:  GOTO  20160 
201'tO  Q(2)  =  Q(2)  -  1 

20143  DEF  SEG  =  SHIAOO:  CALL  RNGENIU!) 

20145  S2INC!  =  FNEXPONT!(U! ,  U2!  )     'compute  new  server  2  service  duration 
20150  NSCT!t2)  =  S2INC!  +  CTIME!      'compute  new  service  completion  time 
20160  GOSUB  26000    "generate  new  exit(2)  indicator  variable 

20170  GOTO  5050     'return  to  clock  module 

20171  ' 

22000  '================  SERVER  2  PLOTTING  SBRT  /  SC2M  ======================= 

22010  LINE  (276,36)  -  ( 290 ,50  )  ,0,BF   'erase  circle  3  server  2 
22020  IF  EXIT! 2)  =  1  THEN  GOSUB  23000:  GOTO  22060   "flash  arrow  4 
22030  INDEX  =  VAL( TRANSFER2$ )  +1   'increment  feedback  counter 
22035  GOSUB  23100   "flash  arrow  5 

22040  IF  STAT(l)  =  0  THEN  X=140  :  Y=36  :  GOSUB  27000  :  GOTO  22055 
22050  IF  Q(l)  <  =7  THEN  X  =  116  -  Q(l )  »  16  :  Y=28  :  GOSUB  27000 
22055  GOSUB  28000:  PUT(15,76),  ARROWS,  PSET  'plot  white  arrow  5 
22060  IF  Q(2)  =  0  THEN  RETURN  'return  to  line  20090 

22069  "erase  ball  at  server  1  and  place  at  server  2 

22070  GET(252,28)  -  (266,42),  TEMPCIRC  :  LINE( 252,28  )-( 266,42  )  ,0,BF: 

PUT( 276,36),  TEMPCIRC 
22080  IF  Q(2)  =  1  THEN  RETURN   'return  to  line  20090 
22090  X=236 
22100  IF  Q(2)  <=  6  THEN  K  =  (Q(2)-l)  ELSE  K  =  5 

22109  'shift  all  remaining  customers  one  to  the  right 

22110  FOR  J  =  1  TO  K 

22120     GET(X,28)  -  (X+14,42),  TEMPCIRC 

22130     LINE(X,28)  -  (X+14,42),  0,  BF 

22140     PUT(X+16,28),  TEMPCIRC 

22150     X=X-16 

22160  NEXT 

22170  IF  Q(2)  <=  6  THEN  RETURN 

22179  'go  to  overflow  queue  and  bring  next  customer  into  visible  queue 

22180  LO0PC0UNT2$  =  MID$(STRG2$,  ( LEN( STRG2$  )-6  )  ,1 ) 
22190  X  =  172  :  Y  =  28 

22195  INDEX  =  VAL( L00PC0UNT2$ )  +  1 

22200  GOSUB  27000     'bring  circle  in  from  queue  overflow 

22205  RETURN    'return  to  line  20090 

22206  ' 

23000  •================   FLASH  ARROW  4  SBRT  /  SC2M  ========================== 

23010  PUT( 295,43),  ARR0W4C,  PSET   'flash  purple  arrow  4 

23015  SOUND  3000,  DUR2 

23020  FOR  I  =  1  TO  1000 

23030  NEXT 

23040  PUT  (295,43),  ARR0W4,  PSET  'plot  white  arrow  4 
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25050 
23051 
23100 
23110 
23120 
23130 
23135 
231<+0 
231^1 
23500 
23501 
23510 
23520 

23530 
2'+990 
24991 
2'+992 
2<^993 
2499"* 
25000 
25004 
25005 
25010 
25015 
25016 
26000 
26004 
26005 
26010 
26015 
26016 
27000 
27010 

27020 
27030 
27040 
27050 
27060 
27070 
27080 
27090 
27100 
27110 
27111 
28000 
28010 
28020 
28030 
28031 
29000 
29005 
29010 
29020 
29021 
30000 
30008 
30009 
30010 
30015 
30020 
30025 
30030 
30035 
30040 
30045 
30050 
30055 


RETURN 


'return  to  line  22020 


•===  ===:==  =  =  r=====   FLASH  ARROW  5  SBRT  /  SC2M  ====  = 

PUT( 15,76),  ARR0W5C,  PSET    'flash  purple  arrow  5 

FOR  1=1  TO  500 

NEXT 

SOUND  500,DUR1 

RETURN    'return  to  line  22040 


>========  FREQUENCY  OF  SERVICE  BY  SERVER  2  SBRT  /  SC2M  =============== 

'this  sbrt  used  in  calculation  of  S2 

FOR  1=0  TO  9 

IF   VAL(TRANSFER2$)    =   I   THEN  SERVBY2(I+1)    =   SERVBY2(I+1)    +   1: 

SERVBYKI  +  l)    =   SERVBYKI  +  l)    +    1:    RETURN 
NEXT 

^f^^^9(i(9(J(XXMXXKK^}t^WM^Jt'XKXXXJ(XMJ^WXMKMMKitXXi(JfXwXw)iXX^(w9( 

***********  SUBROUTINES  USED  BY  MORE  THAN  ONE  MODULE  *************** 
******************************************************************** 


========  PROB  EXIT  AFTER  SERVER  1  SBRT  /  CM,  SCIM, 

exit  after  service  at  server  1:  l=yes  0=no 

DEF  SEG  =  8H1A00:  CALL  RNGENCU!) 

IF  U!  <=  PI!  THEN  EXIT(l)  =  1  ELSE  EXIT(l)  =  0 

RETURN   'return  to  line  5034  or  15140 

>========  PROB  EXIT  AFTER  SERVER  2  SBRT  /  CM,  SC2M  = 

'exit  after  service  at  server  2:  l=yes  0=no 
DEF  SEG  =  8H1A00:  CALL  RNGENCU!) 
IF  U!<=  P2!  THEN  EXIT(2)  =  1  ELSE  EXIT(2)  =  0 
RETURN   'return  to  line  5050  or  20170 


•===============  QUEUE  OVERFLOW  SBRT  /  SCIM,  SC2M  ================= 

ON  INDEX  GOTO  27020,  27030,  27040,  27050,  27060,  27070,  27080, 

27090,  27100,  27110 
RETURN  'return  to  line  17180  or  22040  or  22205 
RETURN 
RETURN 
RETURN 
RETURN 
RETURN 
RETURN 
RETURN 
RETURN 
RETURN 


PUTtX,Y),  CIRCO 

PUT(X,Y),  CIRCl 

PUT(X,Y),  CIRC2 

PUTtX,Y),  CIRC3 

PUT(X,Y),  CIRC4 

PUT(X,Y),  CIRC5 

PUT(X,Y),  CIRC6 

PUT(X,Y),  CIRC7 

PUT(X,Y),  CIRCS 

PUTIX.Y),  CIRC9 


DELAY  SBRT  /  SCIM,  SC2M 


FOR  I  = 

NEXT 

RETURN 


1  TO  500 


•============  CAPTURE  KEYBOARD  INPUT  SBRT  /  MPM,  CM 

DEF  SEG:  POKE  8H6A,0: 

RESPONSE$  =  INKEY$:  IF  RESPONSE*  =  ""  THEN  29010 

RETURN 


•=============-  PRINT  STATISTICS  SBRT 

'sbrt  prevents  displaying  number  that 
IF  (*R1!  <  .001  AND  R1!>0)  OR  (Rl!  > 
LOCATE  14, Y:  PRINT  USING  FMT2$>  Rl! 
IF  ((R2!  <  .001  AND  R2!>0)  OR  (R2!  > 
LOCATE  15, Y:  PRINT  USING  FMT2$>  R2! 
IF  (ISl!  <  .001  AND  S1!>0)  OR  (SI!  > 
LOCATE  16, Y:  PRINT  USING  FMT2$>  SI! 
IF  ((S2!  <  .001  AND  S2!>0)  OR  (S2!  > 
LOCATE  17, Y:  PRINT  USING  FMT2$>  S2! 
IF  ((LI!  <  .001  AND  L1!>0)  OR  (LI!  > 
LOCATE  18, Y:  PRINT  USING  FMT2$>  LI! 
IF  ((L2!  <  .001  AND  L2!>0)  OR  (L2!  > 


/  MPM, 
exceed 

CM  ====== 

formats 

._. 

9999!  )) 

THEN 

Rl! 

= 

9999 

999 

9999!  )) 

THEN 

R2! 

= 

9999 

999 

9999!  )) 

THEN 

SI! 

= 

9999 

999 

9999!  )) 

THEN 

S2! 

= 

9999 

999 

9999!  )  ) 

THEN 

LI! 

= 

9999 

999 

9999!  )) 

THEN 

L2! 

= 

9999 

999 
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30060  LOCATE  19, Y:  PRINT  USING  FMT2$>  L2! 

20065  IF  ((L!  <  .001  AND  L!>0)  OR  (L!  >  9999!))  THEN  L!  =  9999.999 

30070  LOCATE  20, Y:  PRINT  USING  FMT2$>  L! 

30075  IF  ((Wl!  <  .001  AND  H1!>0)  OR  (HI!  >  9999!))  THEN  Wl!  =  9999.999 

30080  LOCATE  21, Y:  PRINT  USING  FMT2$>  Wl! 

30085  IF  ((W2!  <  .001  AND  H2!>0)  OR  (W2!  >  9999!))  THEN  H2!  =  9999.999 

30090  LOCATE  22, Y:  PRINT  USING  FMT2$>  W2! 

30095  IF  ((W!  <  .001  AND  W!>0)  OR  (W!  >  9999!))  THEN  H!  =  9999.999 

30100  LOCATE  23, Y:  PRINT  USING  FMT2$>  W!j 

30105  IF  ((Z!  <  .001  AND  Z!>0)  OR  IZ!  >  9999!))  THEN  Z!  =  9999.999 

30110  LOCATE  24, Y:  PRINT  USING  FMT2$>  Z! > 

30130  RETURN   "return  to  line  <t936  or  9210 
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